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Abstract. We present the results of an extensive search for periodic and irregular variable pre-main sequence (PMS) stars 
in the young (2-4 Myr) open cluster NGC 2264, based on photometric monitoring using the Wide Field Imager (WFI) on the 
2.2 m telescope on La Silla (Chile). In total, about 10600 stars with Ic magnitudes between 9.8 mag and 21 mag have been 
monitored in our 34' x 33' field. Time series data were obtained in the I c band in 44 nights between Dec. 2000 and March 2001; 
altogether we obtained 88 data points per star. Using two different time series analysis techniques (Scargle periodogram and 
CLEAN) we found 543 periodic variable stars with periods between 0.2 days and 15 days. Also, 484 irregular variable stars 
were identified using a^ 2 -test. In addition we have carried out nearly simultaneous observations in V, Rc and a narrow-band 
Hor filter. The photometric data enable us to reject background and foreground stars from our sample of variable stars according 
to their location in the I c vs (R c - Ic) colour-magnitude and (R c - Ha) vs (R c - Ic) colour-colour diagrams. We identified 405 
periodic variable and 184 irregular variable PMS stars as cluster members using these two different tests. In addition 35 PMS 
stars for which no significant variability were detected could be identified as members using an Ha emission index criterion. 
This yields a total of 624 PMS stars in NGC 2264, of which only 182 were previously known. Most of the newly found PMS 
stars are fainter than Ic — 15 mag and of late spectral type (SM2). We find that the periodic variables, as a group, have a smaller 
degree of variability and smaller Ha index than the irregular variables. This suggests that the sample of periodic variables 
is biased towards weak-line T Tauri stars (WTTSs) while most of the irregular variables are probably classical T Tauri stars 
(CTTSs). We have quantified this bias and estimated that the expected fraction of WTTSs among PMS stars in the cluster 
is 77%. This is relatively close to the fraction of WTTSs among the periodic variables which is 85%. We also estimated the 
total fraction of variables in the cluster using only two well selected concentrations of PMS stars called NGC 2264 N & S in 
which we can easily estimate the total number of PMS stars. We find that at least 74% of the PMS stars in the cluster with 
Ic < 18.0 mag were found to be variable (either periodic or irregular) by our study. This number shows that our search for PMS 
stars in NGC 2264 through extensive and accurate photometric monitoring is very efficient in detecting most PMS stars down 
to at least Ic = 18.0 mag. 

Key words, open clusters and associations: individual (NGC 2264) - stars: pre-main sequence, rotation, activity, spots - meth- 
ods: data analysis, time series analysis, periodogram 



1. Introduction 

The angular momentum problem is one of the longest stand- 
ing conundrums in astronomy. Briefly stated, molecular cloud 
cores have four to six orders of magnitude too much angular 
momentum to be incorporated into a star. Without shedding 
their excess angular momentum, stars could not form. When 
and how they rid themselves of angular momentum is still an 
open question. It is intimately connected to the issue of disk 
formation and evolution, to the formation of binary systems 
and, during the pre-main sequence (PMS) phase, to the mag- 
netic interaction between stars and disks which likely regu- 
lates accretion and drives outflows and jets (e. g. Bodenheimer, 
fl989l Bouvier et al. [l997l Stassun et al. ll9^9l Mathieu l2003l . 
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Observationally, our best hope of constraining the angular 
momentum evolution of PMS stars is to obtain stellar rotation 
rates for objects at a variety of masses and ages. Fortunately, 
it is possible to do this photometrically because PMS stars 
have strong surface magnetic fields which produce large cool 
spots (e. g. Feigelson & Montmerle 1999 1. It is often the case 
that the spot pattern is sufficiently asymmetric and stable that 
the rotation period can be found by photometric monitoring 
(e. g. Rydgren & Vrba [19831 Herbst et al. Il994l Stassun et 
al. 1999). The advantage of this photometric method compared 
with v sin i measurements is that it yields directly the rotation 
period independent of the inclination of the star and it can be 
very accurate even for slow rotating stars. Herbst et al. ([2001 
2002) have used this method and obtained rotation periods for 
369 T Tauri stars (TTSs) in the ~ 1 Myr old open cluster Orion 
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Nebular cluster (ONC). They found that the period distribu- 
tion for higher mass stars (i. e. M > 0.25Mq) is bimodal with 
peaks at 2 and 8 days. The bimodality is interpreted as an effect 
of disk-star interactions in PMS stars: Slow rotators have been 
magnetically locked to their disks which prevent them spinning 
up with increasing time or equivalently with decreasing radius. 
Magnetic coupling to an accretion disk was first proposed more 
than a decade ago as a dominant braking mechanism for PMS 
stars rCamenzind fT990l Konid fT99"Tl Shu et al. fT994l i. 

However, we also mention a principal limitation of the 
photometric method for studying the rotational periods among 
TTSs. As we will show the photometric method can much more 
easily measure the periods among the weak-line T Tauri stars 
(WTTSs) but only for a fraction of the classical T Tauri stars 
(CTTSs) due to their higher irregular variability. This "photo- 
metric noise" often prevents the detection of a periodic signal 
in the photometric data. Therefore any photometric monitor- 
ing of a young open cluster (like the ONC or NGC 2264) will 
be biased towards WTTSs in the sense that for these stars the 
fraction of measurable periods will be higher. 

From the results of the rotational study in the ONC several 
questions naturally arise: 1) Is the period distribution similar 
in other clusters (i. e. does environment play a role)? 2) How 
many PMS stars interact with their disks and how strongly does 
this affect the angular momentum evolution? 3) How does the 
period distribution evolve with time (i. e. do most stars decou- 
ple from their disks on time scales of the age of an ONC star 
and spin up with conserved angular momentum or do disks and 
disk-locking persist for longer time scales)? To answer these 
questions it is necessary to measure rotation periods of large 
samples of PMS stars in clusters with different ages. 

Aside from the ONC, the open cluster NGC 2264 is per- 
haps the best target for a detailed rotational study, since it is 
sufficiently nearby (760 pc, Sung et al. U997l l. fairly populous, 
and with an estimated age of 2 - 4Myr (Park et al. 2000 1 it is 
about a factor of 2 - 4 older than the ONC. 

Prior to our study a few monitoring programs have been 
carried out in NGC 2264 (e.g. Kearns & Herbst IT9"9"8T l which 
altogether yielded only about 30 published rotation periods for 
PMS stars. Vogel & Kuhi [T98T1 and Soderblom et al. (JT999J 
reported rotational velocities (v sin for in total about 60 low 
mass NGC 2264 stars. An early variability study has been car- 
ried out by Nandy & Pratt (119721 who identified 26 optical 
variable stars in the cluster. Neri, Chavama, & de Lara ( 1993 1 
reported optical and near-infrared photometry for 50 stars in 
NGC 2264 and detected for 70% of these stars optical variabil- 
ity. Due to the limiting sensitivity of these studies most of the 
stars have spectral types earlier than K5. 

About 200 PMS stars and some 20 PMS candidates have 
been identified in NGC 2264 prior to our study using a variety 
of methods including Ha spectroscopy (Herbig 1954 Ogura 
119841 Marcv lT980l Rvdgren fT979> . Ha narrow band photom- 
etry (Sung et al. lT997l Park et al. l2000l or ROSAT X-ray flux 
measurements (e.g. Flaccomio et al. 2000 1. A PMS member- 
ship catalogue that summarises (most) these results is available 
in Park et al. d2000l . 

Vasilevskis, Sanders, & Balz ( 1965) have carried out a 
proper motion study of 245 stars in the region of the clus- 



ter. However, the suggested membership probabilities must be 
handled with caution because several stars with strong Ha- 
emission which clearly must be members have low member- 
ship probability while other (non-emission) stars with high 
membership probability are clearly no members (see e. g. 
the Ha-emission study by Rydgren 1979 1. King (1998 1 and 
Soderblom et al. (1999 1 have in total determined lithium abun- 
dances of 35 cluster members. Recently, circumstellar disk 
candidates have been identified in NGC 2264 by Rebull et al. 
(120021 . 

In this paper we present the results of an extensive photo- 
metric monitoring program in NGC 2264, in which we moni- 
tored ~ 10600 stars over a very broad magnitude range (9.8 < 
Ic ^ 21). This extensive monitoring program allowed us to dis- 
cover about 400 new PMS stars in the cluster and in addition 
we could measure the rotation periods of 405 PMS stars down 
to the substellar limit. 

In the following section we describe the observational and 
data reduction strategy. Sect.[3]describes the methods we have 
used for obtaining absolute and relative photometry. In Sect. |4] 
we describe the time series analysis. In Sect. [5] we discuss in 
detail how we selected the PMS stars among the periodic and 
irregular variables we have found. The nature of the stellar vari- 
ability is examined in Sect. |6] The identification of additional 
PMS members for which no variability could be detected is de- 
scribed in Sect. [71 In Sect.[8]we estimate the fraction of variable 
PMS stars and the completeness level of our PMS sample. In 
Sect.|5]we investigate if our sample of periodic variables is rep- 
resentative for all cluster members. Our major conclusions are 
summarised in Sect. [K)] The obtained extensive rotational pe- 
riod data and possible consequences for the disk-locking sce- 
nario will be discussed in a separate paper (see Lamm et al. 
2003 hereafter Paper II). 

2. Data acquisition & reduction 

The photometric monitoring program described here was car- 
ried out in the Ic band on 44 nights during a period of two 
month between 30 Dec. 2000 and 01 Mar. 2001 with the Wide 
Field Imager (WFI) on the MPG/ESO 2.2 m telescope on La 
Silla (Chile). The WFI consists of a mosaic of four by two 
CCDs with a total array size of 8 k x 8 kk. The field of view 
is 34' x 33' and the scale is 0.238"/pixel. In Fig.[2we show a 
500 sec Rq exposure of our observed field from which the po- 
sitions of the eight CCDs is evident. To avoid highly saturated 
images and light scattering from the very bright star S Mon 
(V — 4.7 mag) this star was located near the northern end of 
the 23'.'3 (96 pixels) wide central gap between two chips. The 
central position on the sky was close to RA(2000) = 6 h 40 m 59 s 
and DEC(2000) = 9° 38' 59" for most frames. The typical de- 
viation from this nominal position is 2". About 5% of the im- 
aged area is lost due to the gaps between the chips and the 
small dithering of the frames. In order to increase the dynam- 
ical range of the observations, three consecutive exposures of 
5 sec, 50 sec and 500 sec were taken with the Ic filter. This se- 
ries of three exposures is defined here as one data point. In total 
we obtained 110 data points and between 1 and 18 data points 
per night. The observing time distribution of our time series is 
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Fig. 1. Reproduction of a 500 sec Rc band exposure of our observed field. The area imaged by the individual chips a - h is quite 
evident. The bright star near the northern end of the gap between chip b and c is S Mon (V=4.7). In the second chip g the famous 
cone nebular is visible, which is hardly evident on our Ic band images as is the case for the nebulosity in the centre of chip e. 
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Fig. 2. Distribution of observation times of the Ic band time 
series (between 30 Dec. 2000 and 01 Mar 2001). The time (ex- 
pressed as a fraction of a day) as a function of the modified 
Julian day (MJD) is shown. 

shown in Fig. [3 The typical seeing (measured by the FWHM 
of the PSF) in these images was of the order of 0(!%-Y!2. 

In addition to the Ic observations we observed the cluster 
through V and Rq filters on six nights during the 2000/2001 
season. In Rc the exposure times were 5 sec, 50 sec and 
500 sec, while the exposure times through the V filter were 
set to 5 sec, 60 sec and 720 sec. Ic observations were obtained 
directly before or after an observation in V or Rc- This allows 
us to determine the (V - Ic) and (Rc - Ic) colours from nearly 
simultaneous measurements in the two filters, i. e. the colours 
should not be affected by variability to any substantial extent. 

In order to improve the absolute photometry, additional V, 
Rc and Ic WFI data were obtained between October 2001 and 
March 2002 using the same exposure times. In order to search 
for Ha emission, three consecutive Ha observations of 150 sec, 
1200 sec and 1200 sec exposure time were taken directly be- 
fore or after the observations in Rc in order to obtain nearly 
simultaneously observations in these two filters for a determi- 
nation of the (Rc - Ha) colour. The used Ha filter has a central 
wavelength A = 6588 A and FWHM = 74 A. A detailed list of 
observations in filters other than Ic is given in Tabled 

The image processing was done for each of the eight 
WFI chips separately using the standard IRAF tasks. A one- 
dimensional bias was subtracted from each frame using the 
overscan region in each frame. A small two-dimensional resid- 
ual remained which was removed using zero integration time 



Table 1. Dates of observations with the V, Rc and Ha filter for 
the two observing seasons 2000-2001 and 2001-2002. Listed 
are the number of data points (i. e. three images of three differ- 
ent exposure times) for each night in a given filter. The obser- 
vations in V and Rc were take directly before or after a set of 
three images in Ic- 



Date of 


data points 


in 


observation 


V 


Rc 


Ha 


9 Jan 2001 


1 






10 Jan 2001 


1 






11 Jan 2001 


1 






12 Jan 2001 


1 


1 




15 Jan 2001 




1 




16 Jan 2001 




1 




25 Nov 2001 


1 


1 


1 


11 Dec 2001 


1 






4 Jan 2002 


1 


1 


1 


13 Jan 2002 


1 






17 Jan 2002 


1 






2 Mar 2002 


1 


1 


1 


total 


10 


6 


3 



(bias) frames. For images through V, Rc and Ha filters the vari- 
able sensitivity across the chips was corrected using the median 
combination of typically 10- 15 twilight flats taken in the be- 
ginning and the end of two or three different nights. In Ic we 
used illumination-corrected dome flats for flatfielding the sci- 
ence images. This was necessary because the Ic twilight flats 
showed interference fringes caused by narrow band emission 
from the Earth's atmosphere. The illumination-corrected flats 
were created in the way described by Bailer- Jones & Mundt 
(2001 1. Since the science images in the /c-band were affected 
by fringing too, it was essential that these fringes be removed. 
We note that fringing is an additive effect and therefore the 
fringe pattern has to be subtracted from the science images. 
Since the fringe pattern was found to be stable over a given 
observation run we created only one fringe image for each 
chip per observing season (for more details see Bailer-Jones 
& Mundt. l200TV 

3. Photometry and Astrometry 

Our goal was to perform relative and absolute photometry for 
every star in our field with a sufficient signal-to-noise ratio, i. e. 
S IN - 10 or better. The relative and absolute photometry is 
described in the following subsections separately. As a com- 
mon first step we had to identify the positions of all sources 
in our field. Therefore we first created a source list for all 8 
chips employing the DAOFIND task, with a separate list for 
the 5 sec, 50 sec and 500 sec exposures. For this source search 
we selected for each of the three 5 sec, 50 sec and 500 sec ex- 
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Table 2. Number of detected sources in the various chips of the 
reference images. 



Exposure 




Number of sources in 


Chip 








time 


a 


b c d e 


f 


g 


h 


total 


5 sec 


505 


268 445 311 405 


278 


286 


334 


2832 


50 sec 


864 


447 586 883 814 


546 


432 


648 


5220 


500 sec 


2220 


900 1187 1535 1754 


1054 


654 


1645 


10940 



posures one frame with a position very close to the nominal 
central location of all frames. These images we designate here 
as reference frames. The number of sources detected in each 
chip of these reference images is listed in table |5] The coor- 
dinates of each source in all other frames were calculated by 
applying a derived linear offset between the reference images 
and these frames. The offsets were calculated for each chip sep- 
arately by measuring the pixel position of a bright star close to 
the chip's centre in all frames. For the whole observed field the 
offsets were calculated by using eight bright stars. During the 
aperture photometry all sources were re-centered (see below) 
and therefore the offsets calculated in that way are sufficient 
for identifying the stars in the all other frames. 

The IRAF tasks CCXYMATCH and CCMAP were used to 
cross-identify our sources in the reference frames with objects 
in the USNO-A2.0 catalogue (Monet et al. 1998 1 and to calcu- 
late the plate transformations to the sky positions. We typically 
used 70 reference stars (minimum 30, maximum 133) for the 
calculation of the plate solutions in each chip. Using these plate 
solutions the final J2000 coordinates of our sources were cal- 
culated with the task CCTRAN. These coordinates are listed 
in Table 0] The RMS of the deviations between the fit and the 
coordinates of the reference stars is typically 071. 

3.1. Relative photometry 

All the variability studies described here are based on differ- 
ential photometry relative to a set of non-variable reference 
stars. Before the final analysis and after various tests, we re- 
jected some frames from the further analysis. These were usu- 
ally frames that were taken under poor transparency conditions, 
poor seeing or with a bright background due to scattered moon 
light. In total we rejected 22 of the 110 data points before the 
final analysis, so that relative photometry was performed on 88 
remaining data points. 

The DAOPHOT/APPHOT task was used to measure the 
brightness of each object in the source lists. The aperture was 
chosen to be 8 pixels (179) in diameter for all measurements in 
order to maximise the signal-to-noise ratio. The sky was cal- 
culated for each source separately as the median of an annulus 
with an inner diameter of 30 pixel (771) and a width of 8 pixel 
centred on the source. For the calculation of the sky value we 
used a sigma clipping rejection criterion with a 3cr-threshold. 
During the measurement the sources were re-centred. Sources 
with any pixel entering the < 1 % nonlinearity region were re- 
jected. 

The differential magnitudes for all sources were formed as 
follows. First a set of non-variable comparison stars was se- 



Fig. 3. The mean error 6m K \ in a single measurement for the 
three different exposure times of 5s (top), 50s (mid) and 500s 
(bottom) as a function of Ic magnitude. The shaded regions 
indicate which exposure time we used for the further analysis 
of a star with a given magnitude (see text). Note the different 
magnitude scales in the three plots. 

lected from all sources in each chip for each of the three dif- 
ferent exposures. Comparison stars were chosen according to 
the following criteria: 1) present on every of the 88 frames, 2) 
isolated from other sources and the corners of the field (see be- 
low), 3) the photometric error given by APPHOT (i. e. Poisson 
error in the source and scatter in the background) is on average 
less than 0.01 mag. Only stars that passed all of these three tests 
were selected. Possible variables among the candidate compar- 
ison stars were identified by comparing the flux of each can- 
didate with the mean flux of all other candidates at the differ- 
ent epochs. Stars that show the largest variability were rejected 
from the candidate sample and the procedure was performed 
again with the remaining stars. After a few iterations typically 
ten non-variable comparison stars were identified in each chip 
so that the mean standard deviation in the relative light curves 
of the comparisons stars (relative to the other comparison stars) 
was typically <x = 0.009 mag. In the best cases we achieved a 
mean standard deviation of cr = 0.006 mag which is therefore 
the maximum precision we can expect for the relative photom- 
etry. This limit is most likely set by flat fielding errors and not 
by photon noise. 

The selected comparison stars were used for the calcula- 
tion of the relative magnitude m re i(f,) of each source in the field 
at the different epochs f,. The relative magnitude and its error 
Sm K \(ti) were calculated as described in Bailer- Jones & Mundt 
(2001 1. We note that the mean was subtracted from each light 
curve so that i*ke\(ti) = 0. 

In Fig.|3]we show the mean error 5m K \ in a single measure- 
ment as a function of magnitude for the stars measured with 
the three different exposure times (for the magnitude determi- 
nation see section l3~2l . As an example we show the light curves 
of four stars in Fig.0] 

3.2. Absolute photometry 

In order to obtain additional constraints on the nature of the 
observed stars, in particular on their PMS membership, we ob- 
tained absolute photometry in V, Rc and Ic- As we will outline 
below all of our photometry is done relative to secondary stan- 
dards in our observed field. The result of the absolute photom- 
etry is reported in Table0] 

Many of the stars in the field are expected (and found, see 
Section^} to be variable and the peak-to-peak variations of the 
variable stars are in most cases £ 0.2 mag in Ic (see Fig. ll 31 . 
Therefore a single measurement of the magnitude of a star de- 
pends on the phase and amplitude of the star's brightness mod- 
ulation. On the other hand single measurements of the colour 
of a PMS star with cool (or hot) surface spots also differ at dif- 
ferent epochs. That is because the spectral energy distribution 
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Fig. 4. Examples for relative light curves (mean subtracted) of four stars and the resulting phased light curve for each of the four 
stars. The light curves were phased with the period that we found by the Scargle periodogram analysis technique. For each star 
we also list the Number (no.), the mean Ic magnitude and the period P. 

Table 3. Comparison of our photometry with those of other authors. We list the offsets in the sense results of other author minus 
our results. We also show the number of stars which we used for the comparison with other studies. These numbers differ for a 
single publication because not all magnitudes or colours were available for all stars. Stars with close-by neighbours were not not 
used for the calculation of the offsets, since these stars are not separated in the other studies. 



Publication 


No. of stars 


Ic 


No. of stars 


(Rc 


-Ic) 


No. of stars 


(V-Ic) 


Rebull et al.(2002) 


1344 


0.003 + 0.002 


1047 


-0.019 


± 0.002 


767 


-0.010 + 0.002 


Park et al. (2000) 


147 


-0.009 ± 0.012 








147 


-0.041 ± 0.002 


Flaccomio et al. (1999) 


236 


0.070 ± 0.003 


236 


-0.072 


± 0.008 


221 


-0.147 + 0.007 


Sung etal. (1997) 


117 


-0.000 ± 0.002 








116 


0.001 ±0.010 



of the light coming from the area of the spot differs from the 
spectral energy distribution of the light radiated from the other 
part of the stellar atmosphere because of the different effec- 
tive temperatures. The difference between the two energy dis- 
tributions is wavelength dependent and larger for shorter wave- 
lengths. Therefore, if cool spots cause the variability of a star 
the peak-to-peak variation in its light curve increases to shorter 
wavelength, i. e. it is larger in the V-band and smaller in the 
Ic band. Observations of T Tauri stars confirm this expected 
colour behaviour of the stars (e.g. Vrba et al. [1989 Herbst et 
al., [T994i. The relations between the peak-to-peak variations 
AV, ARc and AIq in the different filters are given by the colour 
slopes 

d(ARc-Mc) , _ d(AVc-Mc) 

Sr = and Sv = ■ 

d(AI c ) d(AIc) 

Herbst et al. (1994i calculated these slopes for several WTTSs 
with a mean value of S R = 0.31 ± 0.23 and S v = 0.55 + 0.36. 
Since the typical peak-to-peak variations we found in our target 
stars are of the order of Ale - 0.2 we expect maximal colour 
changes of about A(R C - Ic) - 0.06 and A(V - Ic) - 0.1 but in 
some extreme cases the variations in the colour may be higher. 

To take into account these magnitude and colour changes 
we determined the Ic magnitude and the (V - Ic) and (Rc - Ic) 
colours at different epochs. The colours and magnitudes we list 
in Tableware average values of these different measurements. 
The dates of the employed measurements are listed in Tabled 

For the determination of the average Ic magnitude of each 
star we used only images which were taken during the first ob- 
serving season in Jan. 2001. In this season we got nearly si- 
multaneous measurements in Ic and Rc (i. e. three images in Ic 
with different exposure times followed by three images of dif- 
ferent exposure times in Rc) or nearly simultaneous measure- 
ments in Ic and V at seven different epochs (see Table[0. The 
final averaged (Rc - Ic) colours were calculated from nearly 
simultaneous measurements in Rc and Ic at the six different 
epochs during both observing seasons listed Tabled The aver- 
aged (V— Ic) colours were calculated from nearly simultaneous 
measurements in V and Ic at the 10 different epochs in both 
observing seasons listed in Tabled The transformation of our 
instrumental Vmstr, ^c.instr and /c.instr magnitudes into the true 



(V — Ic) and (Rc - Ic) colours and Ic magnitude is outlined be- 
low. It was done before the averaging and for each epoch, each 
chip and each exposure time separately using a linear transfor- 
mation. 

Since we did not observe any flux standards during any of 
the observing seasons we had to use secondary standard stars 
located in our field for calibrating our measurements. The pho- 
tometric calibration coefficients were determined using mag- 
nitudes and colours of stars in NGC 2264 measured by Rebull 
et al. (private communication). Their photometric data include 
the photometry of PMS stars in the cluster from Rebull et al. 
(2002) and in addition unpublished photometry of foreground 
and background objects in their observed field. This extended 
dataset including the unpublished data, compared to a dataset 
consisting only of PMS stars, has the advantage of a smaller 
fraction of variable stars. Our objects were cross-identified with 
the 2924 objects in the extended Rebull et al. catalogue that are 
located in our field. The 2227 (76%) stars which we could iden- 
tify we used as secondary photometric standard stars. 

We transformed our instrumental magnitudes 7c,insti into the 
Cousins / system by applying the linear transformation equa- 
tions of the form 

I C = a i + b\ X (flc.instr ~ Vinstr) + fc,instr and 
Ic — ai + b2 X (Vinstr - /c.instr) + /c.instr 

to the measurements in Ic which were nearly simultaneous with 
one of the other two filters V or Rc- The coefficients a\ and b\ 
were calculated simultaneously using the secondary photomet- 
ric standard stars by applying a linear least squares fit to the 
points in the (/ c ,Rebuii - ^c,fastr) vs (/? c ,instr - /c.instr) plane, where 
/c.Rebuii is the magnitude in our secondary standard catalogue. 
The coefficients 02 and Z?2 were calculated in the same way us- 
ing the (/c.Rebuii - /c.instr) vs (V instr - / c ,in S tr) plane. 

We typically used between 100 and 200 stars per chip and 
exposure time (minimum 42, maximum 392) for calculating 
the linear transformation coefficients. Thus even though many 
of the standard stars are variable this large number ensures a 
robust transformation. Obvious outliers were rejected before 
the calculation of the coefficients was done. Since the variable 
stars were in different phases of their brightness modulation 
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for a given epoch the related data points scatter around the true 
transformation function. The uncertainties of the transforma- 
tion coefficients results from this scatter. The median values of 
the slopes are b\ = 0.18 and b 2 = 0.10. The typical uncer- 
tainties of the transformation coefficients are 6a i = 0.005 and 
6b i =0.010 for the first and 6a 2 = 0.015 and 6b 2 = 0.010 for 
the second equation. With this method we derived values for Ic 
at 7 different epochs between 9 Jan 2001 and 16 Jan 2001. 

The final Ic magnitude for each star was calculated as the 
average of these 7 different Ic measurements. To improve the 
results we make use of the fact that we know the phase and am- 
plitude of each star at these epochs from the relative photomet- 
ric analysis described in Sect. 13.11 Before the final averaging 
we subtracted the relative magnitude (m re i) of the correspond- 
ing data point in the (mean subtracted) relative light curve from 
the calculated Ic magnitude at each epoch. The mean of the re- 
sulting values is listed in Tabled 

The errors for this final Ic magnitudes were estimated in 
two different ways. First, we calculated the expected Icr error 
<5^c.exp from the measured IRAF errors 1 and the errors of the 
transformation coefficients. Second, we calculated the standard 
deviation <x 7 of the 7 independent Ic measurements. The final 
assigned error to the magnitude is the maximum of this two 
error estimations for each star: 61c = max(6Ic,exp, cr 7 ). 

The colours in the Johnson V and Cousins R, I system were 
calculated in a similar way. The transformation from instru- 
mental colour to the true colour was done using the transfor- 
mation equations 

(R C - Ic) =C!+d\X (7?C,instr ~ /c.instr) Or 
(V - I C ) = Ci+dlX (Vjnstr - /c,insti)- 

Again, the transformation coefficients were calculated using a 
least square linear fit. The median values of the slopes are d\ = 
-0.26 and di = -0.17. With this method we derived ten values 
for (V - Ic) and six for (Rc - Ic)- The final colours for each star 
were obtained by calculating the mean of these values for each 
colour. The errors were derived in the same way as described 
above. 

3.3. Final photometry database 

The absolute photometric calibration described in the previ- 
ous section was done for each chip and exposure time sepa- 
rately. The objects that were measured on the same chip with 
two different exposure times were identified in order make 
sure that they do not appear twice in the final catalogue. For 
stars with Ic > 16.0 mag the colours and magnitudes listed 
in Table 0] were taken from the 500 sec exposures. The mea- 
surements for stars with 12.5 mag < Ic < 16.0 and with 
10 mag < Ic < 12.5 mag were taken from the 50 esc and 5 sec 
exposures, respectively. Thus we list the measurements with 
the highest signal-to-noise for each star and avoid pixel satura- 
tion in images with the best seeing. The magnitude ranges used 

1 The IRAF errors were corrected to somewhat higher values as 
described in section l4~2l 



Table 5. Limiting magnitudes for a signal-to-noise ratio of 20 
and 50 of observations in each of the filters V, Rc and Ic- Also 
listed is the estimated limit for completeness in each filter. 



Filter 


Limiting 


; Magnitude 




20o- 


50cr 


V 


22.1 


21.1 


R c 


21.5 


20.4 


Ic 


20.7 


19.6 



for our data analysis of the three different exposure times are 
also indicted in Fig. [3] 

We compared the results of our photometry with those of 
other authors. In Table[3]we report the mean offsets in the pho- 
tometry between our study and those of other studies. We did 
not find any significant difference between our results and those 
reported by Rebull et al. d2002t . Park et al. d2000t . and Sung 
et al. dl997> . The small differences between the studies can 
be explained by intrinsic stellar variability. However, there is 
a significant offset in the photometry compared the the study of 
Flaccomio et al. (2000). The reason therefor could be different 
transmission curves of the used filters and a different average 
colour of the employed flux standards. 

A small part of the final photometric catalogue is shown in 
Tableland the complete table that contains all 10554 stars is 
available electronically. In Table |4] we also list the (Rc - Ha) 
colour of the stars. Since no Ha standard star measurements 
were available only the instrumental colour is given. The limit- 
ing magnitudes in the different pass-bands for different signal- 
to-noise levels is given in Table[5] 

4. Time series analysis 

Our aim was to check all 10554 stars in the field for both pe- 
riodic and irregular variability. For detecting these two types 
of variability we have used different techniques which are de- 
scribed in the following subsections. 

4.1. Periodic variables 

We used two periodogram analysis techniques to search for sig- 
nificant periodicity in the light curves of each of the 10554 
monitored stars. Those stars which have fewer than 15 data 
points in their light curves were rejected from this analysis. 
96% of the remaining 10503 stars have 80 or more data points 
in their light curves. In the following we describe the two used 
periodogram algorithms separately. 

4.1 .1 . The Scargle Periodogram 

First we used the periodogram technique for unevenly sampled 
data described by Scargle ( 1982 1 and Home & Baliunas ( 1986 1 
to search for significant periodicity in all monitored stars. The 
algorithm calculates the normalised power Pn(co) for a given 
angular frequency co = 2nv and identifies the location of the 
highest peak in the calculated periodogram of each star. Finally, 
the light curve of each star is phased to the period according to 
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Table 4. Photometric data, the complete table available electronically. The column called "vari" lists the results of the time series analysis: pv 
to irregular variables (see also the comments in the text), cr is the standard deviation of the light curves from which outliers were rejected (see ti 
peak-to-peak variation in the light curves. 



Star a(J2000) (5(72000) 



V-Ir 



Re - Ic 



Rc - Ha err 



ptp SpT name cross id 



6000 
6001 
6002 
6003 
6004 
6005 
6007 
6008 
6010 
6011 
6012 
6013 
6014 
6015 
6016 
6018 
6019 
6020 
6021 
6022 
6023 
6024 
6025 
6026 
6027 
6028 
6029 
6030 



6:41:8.16 
6:41:8.19 
6:41:8.21 
6:41:8.21 
6:41:8.27 
6:41:8.28 
6:41:8.38 
6:41:8.43 
6:41:8.52 
6:41:8.53 
6:41:8.55 
6:41:8.56 
6:41:8.57 
6:41:8.58 
6:41:8.63 
6:41:8.65 
6:41:8.65 
6:41:8.70 
6:41:8.80 
6:41:8.81 
6:41:8.84 
6:41:8.84 
6:41:8.85 
6:41:8.85 
6:41:8.86 
6:41:8.89 
6:41:8.90 
6:41:8.92 



9:30:4.4 
9:36:56.8 
9:34:9.9 
9:38:30.9 
9:30:23.4 
9:38:14.7 
9:40:36.7 
9:37:53.0 
9:32:51.5 
9:23:6.6 
9:31:6.5 
9:42:52.1 
9:30:4.2 
9:28:14.3 
9:53:47.3 
9:39:11.9 
9:40:14.1 
9:53:25.3 
9:27:54.0 
9:23:43.9 
9:48:48.6 
9:53:1.8 
9:45:28.0 
9:52:48.3 
9:46:1.6 
9:29:45.7 
9:29:29.8 
9:29:22.6 



18.37 
16.91 
15.16 
16.56 
15.56 
13.10 
16.10 
16.70 
19.13 
20.54 
18.02 
14.57 
17.20 
19.94 
17.74 
20.31 
16.67 
20.28 
16.92 
12.91 
21.58 
15.38 
18.83 
20.57 
12.31 
19.50 
19.99 
19.59 



0.02 
0.01 
0.01 
0.01 
0.02 
0.01 
0.01 
0.01 
0.02 
0.06 
0.03 
0.02 
0.01 
0.06 
0.01 
0.08 
0.01 
0.06 
0.04 
0.02 
0.19 
0.03 
0.08 
0.18 
0.02 
0.11 
0.08 
0.06 



4.06 
3.32 
3.83 
2.45 
0.73 
2.64 
3.45 
2.34 
0.90 
4.00 
1.91 
3.35 
-0.50 
3.06 



-0.43 
1.83 



0.71 
1.28 
2.16 
2.45 



0.06 
0.01 
0.08 
0.09 
0.01 
0.01 
0.03 
0.03 
0.04 
0.14 
0.01 
0.02 
0.04 
0.03 



3.06 0.01 



0.38 
0.02 



2.44 0.02 



0.01 
0.05 
0.09 
0.14 



2.31 
2.06 
1.69 
1.93 
1.40 
0.35 
1.48 
2.03 
1.10 
0.44 
2.42 
1.01 
1.97 
0.00 
1.95 
2.29 
1.76 

1.35 
0.98 

1.34 
2.90 

0.38 
0.71 
0.82 
1.12 



0.02 
0.01 
0.01 
0.01 
0.03 
0.01 
0.01 
0.01 
0.01 
0.04 
0.04 
0.01 
0.02 
0.05 
0.01 
0.15 
0.01 

0.03 
0.02 

0.02 
0.04 

0.01 
0.07 
0.04 
0.05 



-2.93 
-3.32 
-3.08 
-2.88 
-2.29 
-3.38 
-3.13 
-2.95 
-2.54 

-3.36 
-3.23 
-2.74 

-2.95 

-3.12 

-3.01 
-3.13 

-2.81 
-2.02 



0.07 
0.02 
0.01 
0.04 
0.08 
0.01 
0.01 
0.13 



0.16 
0.01 
0.01 



0.04 
0.02 

0.03 
0.82 



-3.36 0.02 



IV 

pv 

iv 
iv 

pv? 
pv 



0.02 pv?/iv 



pv 
pv 
pv 



0.02 pv? 
0.01 pv 



pv 



pv 



0.073 
0.053 
0.016 
0.074 
0.105 
0.010 
0.015 
0.015 
0.154 
0.126 
0.211 
0.062 
0.016 
0.182 
0.013 
0.056 
0.010 
0.068 
0.058 
0.073 
0.331 
0.217 
0.021 
1.107 
0.014 
0.238 
0.108 
0.151 



0.26 
0.17 
0.05 
0.27 
0.35 
0.04 
0.06 
0.05 
0.55 
0.52 
0.72 
0.18 
0.05 
0.70 
0.05 
0.21 
0.04 
0.25 
0.24 
0.26 
1.18 
0.72 
0.07 
4.26 
0.05 
1.07 
0.46 
0.70 



G5V 



K7V 



K4V 1 



G5V 1 



R3588 
R3591 
R3590 

Y3607, R3592, MX: 

R3594 

R3595 



Y3618, R3598, S3i: 

R3600 
R3599 

Y3650, R3602, P14' 
R3603 

Y3628, R3606, MX: 



the frequency of the highest peak. In order to decide whether 
there is a significant periodic signal of this period in the light 
curve, the height of this peak has to be related with a false alarm 
probability (FAP) which is the probability that a peak of this 
height is due simply to statistical variations, i. e. noise. 

The standard procedure for determining the FAP is to use 
simulated light curves created with Monte Carlo simulations. It 
is necessary that the time sampling of the simulated light curves 
is identical with that of the measured light curves. The sim- 
ulated light curves represent non-variable stars and therefore 
contain only a intrinsic scatter in the photometry, i. e. noise. For 
each simulated light curve the periodogram has to be calculated 
and the power of the highest peak in that periodogram has to 
be determined. After various simulations (typically ~ 10000) 
the cumulative distribution of the power of the highest peak is 
used to determine the FAP: The FAP of a given power is set 
to be the fraction of simulated non-variable stars which have 
a highest peak power that exceeds P^. If we have simulated 
for example 10000 light curves, the 1% FAP power P PAP=i% 
is the power which was exceeded by the highest peak in 100 
simulations. 

The critical point in this procedure is the simulation of non- 
variable or at least non-periodic stars. The simplest approach to 
doing this is to assume that the data points are statistically in- 
dependent of each other. This is only true if the typical time 



scale for intrinsic variations of the sources is not larger than the 
(typical) time difference of the data points. Strictly speaking 
this assumption is not valid for our sample since the signifi- 
cant time scale for variations in PMS stars is ~ 1 d (e. g. Herbst 
et al. 1 19941 1 and in some nights our data points were obtained 
at closely spaced intervals of much shorter length (~ 0.5 h). 
Therefore we calculated the FAP using different types of simu- 
lations. 

First we make the simplified assumption of uncorrelated 
data points and created light curves of pure noise using Monte 
Carlo simulations. The time sampling was chosen to be the 
same as in the real light curves. Gaussian distributed random 
variables with zero mean and standard deviation cr were as- 
signed to these dates. We set cr = 0.01 mag since the mean 
standard deviation in the light curves of the (non-variable) com- 
parison stars we used for calculating the relative light curves is 
of the same order. Therefore we expect a standard deviation 
of cr = 0.01 in the light curve of a non-variable star. After 
10000 simulations we calculated a 1% FAP power value of 
^fap=i% = 10.2 from the cumulative distribution of the highest 
peak as described above. 

To allow for correlations between the data points a second 
type of simulations was carried out. Instead of simulating pure 
Gaussian noise, the measured light curves of the monitored 
stars were used. For the simulations we kept the time sampling 
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the same as in the measured light curves. The relative mag- 
nitudes were randomly reassigned to the real dates. From the 
cumulative distribution of the maximum power we calculated a 
1%FAP power value of f fap=i% — 9.8. 

To be more conservative the highest power value of the two 
FAP simulations was used to define a cutoff level for the de- 
tection of periodic variables. Out of the 10554 analysed stars 
1 192 were brighter than Iq < 19.5mag and had a peak power 
Pn > 10.2. As an example we show in Fig. |4]the complete 
light curve and the resulting phased light curve we found for 
four stars. In Fig.|5]we show a sample of phased light curves at 
different power levels. 

4.1.2. The CLEAN Periodogram 

The Scargle periodogram technique makes no attempt to ac- 
count for the observational window function W(v), i. e. some 
of the peaks in the Scargle periodograms are normally a re- 
sult of the data sampling. This effect is called aliasing and 
even the highest peak could be an artifact. The CLEAN pe- 
riodogram technique by Roberts et al. ( 1987 1 tries to overcome 
this shortcoming of the Scargle periodogram. The observed 
power spectrum |Z)(v)| 2 is given by the convolution of the (com- 
plex) Fourier transformation F(v) of the continuous signal from 
the star and the spectral window function W(v) of the data set: 

D(y)=F(y)e>W(v) (1) 

(for a detailed description see Roberts et al. 119871 or Bailer- 
Jones & MundtEOHD- 

Even if the spectral window function is known it is not 
possible to de-convolve D(v) directly. Therefore Roberts et al. 
( 1987 1 modified the CLEAN algorithm that is known from the 
reconstruction of two-dimensional images from interferomet- 
ric data. Based on their publicly available Fortran code we have 
created a C realization of the CLEAN algorithm. A description 
of this algorithm is given in Appendix [A] We calculated the 
CLEAN periodogram for each of the 10554 stars. After we had 
identified the highest peak in each of the periodograms the light 
curves of all stars were phased with the corresponding period 

P CLEAN- 

4.1 .3. The error in the measured periods 

The uncertainty in a measured period is set by two fundamental 
limits. First, it is related to the finite frequency resolution 6v of 
the power spectrum which makes it impossible to distinguish 
between closely separated periods and the uncertainty in a pe- 
riod is dv/2. For a discrete data set the resolution 5v is given by 
the width of the main peak of the window function W(v). If the 
time sampling is not too non-uniform this is related to the total 
time span T of the observations with 5v - l/T (Roberts et al. 
1987 1. Since for small errors 6v/6P dv/dp = 1/P 2 the error 
in the period is given by 



For larger errors (i. e. larger periods) this is only a lower limit 
on the period error. Second, for very short periods when the 



finite time span between the data points is comparable to the 
period the uncertainty in a period is given by the sampling error 
which is related to the typical spacing between the data points. 
Both of these fundamental limits leave their footprints in the 
main peak of the window function. 

Since our time sampling is not very uniform (see Fig. [2) 
we determined the error of the periods directly from the width 
of this peak and for each star separately. The full width at half 
maximum of the main peak is typically (i. e. for a light curve 
with 88 data points) vfwhm - 0.02601/day whichis 1.63 times 
larger than what we would expect from the total duration of the 
observations. This can be explained by the non-uniform sam- 
pling. For a light curve with 88 data points we finally assigned 
an error to the period P which is given by 6P = 0.026 x P 2 . For 
light curves with less data points the factor differs depending 
on the width of the main peak of the window function. 

4.1.4. Final period determination 

For all of the 1 192 stars that have a normalised peak power in 
the Scargle periodogram of P N > 10.2 (i. e. FAP < 1%) both 
the light curve and the light curves phased with the periods 
found by both periodogram techniques were checked by eye 
and obvious wrong detections (e. g. due to light contamination 
from a close-by neighbour) were rejected. 

The highest frequency which may be recovered from a sam- 
ple with data points spaced at intervals of A is the Nyquist fre- 
quency v N = 1 /(2A) (Roberts et al. H987i . Using A = 0. 107 day 
which is the median time difference of the data points in a light 
curve we therefore set a lower limit on the periods and accepted 
only periods with P > 1/vn = 2 A = 0.21 day. 

In total we found 543 periodic variable stars. These stars 
are marked in Table|4]with a pv in the column "van". For 520 
(95.8%) of these periodic variables the periods determined with 
the two different periodogram techniques agree within the es- 
timated errors. In all cases where the periods determined with 
both methods did not agree the two periods were beat period of 
each other. For those stars we assigned one of the two periods 
to the star. We always choose the period for which the scatter 
in the phased light curve was lowest. For 17 of these 23 stars 
we assigned the CLEAN period but we kept the Scargle Period 
in six cases. 

In addition all 543 periodic variables we found are listed in 
a separate Table [6] In this table we give the periods we found 
with both the Scargle and the CLEAN periodogram techniques, 
the period we finally assigned to each star, and the estimated er- 
rors for this period. We also list the periods measured by Kearns 
& Herbst ( 119951 . if available. Out of the 543 stars 501 (91.9%) 
have a peak power in the Scargle periodogram of fV ^ 12.40 
which corresponds to a FAP of 0.1%. 

The locations of the periodic variables in our observed field 
are shown in Fig. Their distribution show that there are 
two concentrations of young stars in NGC 2264 which we call 
NGC 2264 N (north) and NGC 2264 S (south). These concen- 
trations are already known in the literature: Sagar et al. 1988 
report these two points of maximum stellar density of cluster 
members. This clustering agrees with the observed distribu- 
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Fig. 5. Examples for phased light curves using the period Pscargie found by the Scargle periodogram technique. The peak power 
Pn in the periodograms of these examples decreases from the upper left to the lower right panel. The examples in the two top 
rows represent the highest power values (fV - 40, . . . , 42) we found, the examples in the two middle rows represent median 
power values (fV =t 28) and the lower two rows represent the lowest power values we accepted (fV — 11, ... , 14). The peak 
power Pn, the Period Pscargie an d the star's identification number is given for each example. The error bar in the upper right 
corner of each panel indicates the mean photometric error in the light curves. 



tion of molecular gas, Ha-emission stars and early type stars 
in the cluster (Mathieu 1986 and references therein). The stars 
located in these two concentrations will be used in Sect. 15. ll to 
determine the region in the colour-magnitude diagram where 
the PMS stars are located. 

Out of the stars which we rejected after our visual inspec- 
tion another 136 were classified as possible periodic variables. 
These stars are marked with a pv? in the column "van" of 
Table[4]but are not listed in Table[6] because they are not used 
for any further analysis (see Paper II). 



4.2. Irregular variables 

Since many of the PMS stars in our study show non-periodic 
brightness modulations with a peak-to-peak variation up to 
1 mag we defined any variable star for which no period could 
be detected as an irregular variable. 

We searched for this irregular variability in the light curves 
of all stars that are not periodically variable in the magnitude 
range of Ic < 19.5 (which corresponds to S /N > 50) that have 
at least 20 data points. Stars with close neighbours were re- 
jected from this sample because our fixed-size aperture could 
lead to seeing-dependent overlapping with a (bright) neigh- 
bour, and hence imitate variability. Since the impact of a faint 
neighbour is negligible we used a rejection criteria that is mag- 
nitude dependent. If two stars were separated by less than 2'.'5 
from each other, both of them were rejected from the analysis 
if the magnitude difference between the two objects was less 
than 2.0 mag. If the magnitude difference of such two stars is 
larger than 2.0 mag only the fainter star was removed and the 
brighter stars was kept. 

In total we analysed 5927 non-periodic stars. 90.6% of 
these stars have more than 70 data points in the light curve. For 
detecting the irregular variables we used ax 2 test, in which we 
calculate the probability that the deviations in the light curve 
are consistent with the photometric errors, i.e the probability 
that the star is not variable. Therefore we evaluate 



X 



Z 



m K \(j) 
, ^W re i(j) 



(3) 



where m K \{ j) is the relative magnitude and 6m le \( j) is the error 
of the j'th data point in the light curve of a star. The probability 
that the light curve of a non-variable star with N data points re- 



sults in a value for chi-square that exceeds the measured value 
X 2 is given by 2 (Press et al. 1992 1 



Q(x 2 \N) = 



rl N-l £\ poo 

I 2 ' 2 ) _ w 



t (N-3V2 e -, dt 



(4) 



The probability P var j that the star is variable is therefore P vai i = 

1 - Q(x 2 \N). 

We note that the^- 2 test is very sensitive to a over- or under- 
estimation of the errors. For example if the error in the single 
measurements is underestimated by a factor of 2 the value for 
X 2 is overestimated by a factor of 4 and we therefore deduce 
an overestimated probability P ran that the star is variable. In 
order to check whether the errors in the light curves are over- 
or underestimated we made a comparison of the scatter <x and 
the average of the errors 5m K \ = (1/N) Yij=i 5m K \ in the light 
curves, where Sm K \ is given by Eq. (4) in Bailer- Jones & Mundt 
(2001 1. In Fig.|6]we show the scatter in the light curve of each 
star measured in the 500 sec exposures. The dotted line is a 
fit M(Ic) to the median of the mean photometric errors of the 
500 sec exposures shown in the lower panel of Fig. [5] The solid 
line represents the fit S{Ic) of the median of the standard de- 
viation cr in the light curves. The median values for both the 
standard deviation cr and the mean error 6m re \ were calculated 
in different Ic magnitude bins. In both cases we fitted an ex- 
ponential function in the fainter regime and a constant in the 
brighter regime. 

It is evident that the estimated errors 6m re \ are systemati- 
cally underestimated for the 500 sec exposures and the same is 
the case for the 5 sec and 50 sec. Therefore we corrected the 
errors measured for each star by calculating the ratio of the two 
functions for each of the three exposure times to give a cor- 
rection function C(Iq) = S(Ic)/M(Ic). The error of each data 
point in the light curves of a given star was corrected using the 
equation 



CWIreUorO') = C(/ c ) 5m Kl ( j) 



(5) 



where we used the mean magnitude for Ic of the star calculated 
as described in Sect. 13.21 By correcting the errors in this way 
we ensure that the relative distribution of the errors in a given 
light curve is conserved. 

We note that we used a conservative fit for the standard de- 
viation in the sense that the standard deviation may be overes- 
timated by our fit (see Fig.|6j. Since the same is the case for the 
fit M(Ic) of the mean errors this effect is minimised (but maybe 
not completely canceled) for the ratio C(Ic) = S(Ic)/M(Ic). 
However, a conservative estimate of the true error reduces the 



2 Since the time series has been mean subtracted the number of de- 
grees of freedom is N-l. 
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Fig. 6. Standard deviation cr of all stars measured with the 
500 sec exposures as a function of magnitude. The solid line 
represents a fit (S (/c)) to median values of cr. The dotted line 
shows the fit (M(Ic)) to the median of the mean photometric 
error Sm le \ displayed in Fig. [3] The medians for both quantities 
were calculated in equally spaced magnitude bins. 

false detections of the x 1 -test. On the other hand the test is less 
sensitive for low amplitude variations. 

From Eq.[3]it is clear that a single outlier in the light curve 
(e. g. caused by a cosmic ray) could produce a very high value 
for the x 2 - Since we are looking only for persistent variability 
we have applied a sigma clipping algorithm to the light curves 
before the^- 2 -test was performed. This algorithm uses the stan- 
dard deviation cr of a given light curve and removes all data 
points from this light curve that are located more than 2.5cr 
above or below zero (note that the mean value was subtracted 
from the light curves). 

After all stars had been analysed we found that many of the 
detected variable stars are located in the corners of our field. 
The light curves of these stars looked very similar and were 
well correlated with the changing seeing. A possible explana- 
tion for this is that the pixel scale of the WFI is different in 
the corners of the field. A similar result was recently found by 
Koch et al. (private communication) for the WFI. This leads 
to variations in the photometry of the stars located in the cor- 
ners because the flux in the aperture of these stars changes with 
the seeing in a different way from the flux in the apertures of 
the comparison stars (the latter are not located in the corners, 
see Sect. 13. 1> . Therefore non-variable stars can mimic variabil- 
ity. In order to avoid this problem resulting from the variable 
pixel scale we performed the relative photometry again in ex- 
actly the same way as described in Sect. B.ll but we used this 
time a significantly larger aperture diameter of 20 pixels (4'.'76) 
instead of 8 pixels. In this way we made sure that much more 
of the star's light falls into the aperture and the measurement 



is much less dependent on the seeing, i. e. the photometric er- 
rors dominate the changes due to seeing variations. The errors 
in the relative magnitudes were corrected in the same way as 
described above. 

In order to decide whether any star in the field is variable 
we used the results of the ^ 2 -tests based on the measurements 
with both aperture sizes. Only the stars that have a probability 
of P V ari ^ 99.9% in both tests were assumed to be irregular 
variable. We used both aperture sizes and not only the mea- 
surements with the bigger aperture because with the large aper- 
ture we have increased the probability that the measurement of 
the stellar brightness is contaminated by the light of a close- 
by neighbour and the star therefore can mimic variability with 
the changing seeing. This effect is minimised for the photom- 
etry with the smaller aperture radius. On the other hand stars 
that do mimic variability because of the variable pixel scale do 
not pass the^ 2 -tests based on the measurements with a bigger 
aperture radius. Hence, a star that passes both tests is not af- 
fected by the variable pixel scale and light contamination from 
a close-by neighbour and therefore the variability is intrinsic to 
the star. 

In total we found 484 irregular variables out of the 5927 
stars we analysed. The irregular variables are listed in column 
"van" of Table |4] as iv. The spatial positions of the irregular 
variables are shown in Fig.0 We remind the reader that irreg- 
ular variables are defined here as variables with P var i > 99.9% 
for which we could not find any significant period in our peri- 
odograms. 

With the x 1 test we are able to detect variability for stars 
brighter than Iq = 16.0 mag if the standard deviation cr in the 
light curves is larger than 0.02 mag. Since we used images with 
500 sec exposure time for stars fainter than Ic < 16.0 mag, 
the sensitivity of the x 1 test is somewhat better for stars with 
16.0 mag < Ic < 17.25 mag. In this magnitude range we are 
able to detect variability for stars with cr > 0.015 mag. For 
stars with Iq = 1 8 mag we are able to detect variability if 
cr > 0.03 mag. 

Only 48.3% of the periodic variables are variable according 
to our^ 2 -test. It is not surprising that not all periodic variables 
are variable according to this test because the periodogram 
analysis is much more sensitive to low amplitude variations 
than the x 2 test (see Sect. |6J. Furthermore we have adopted 
relatively conservative errors for our^ 2 test (see above). This 
decreases the sensitivity to small amplitude variations. 

We note that the different pixel scale in the corners of the 
field has only weak or no effect on the results of our periodic 
variability study. Since the seeing changes randomly the ampli- 
tudes of this random variation contributes to all powers in the 
power spectrum and not only to a single peak. Only ~ 5 - 10 pe- 
riodic variables are located in the affected regions of the field. 
However, the variations in the corners of the field might have 
caused us to miss detecting a small number periodic variables. 

5. PMS Membership 

The periodic and irregular variables we found in our field are 
not necessarily all PMS stars, although variability is probably 
one of the best indicators for youthfulness. Nevertheless there 
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Fig. 7. The locations of all 543 periodic (o) and 484 irregu- 
lar (a) variables from our survey (including non-PMS stars). 
Those variable stars which were previously known as PMS 
stars are indicated by a vertical line ( <|>, 4 respectively) while 
previously known PMS members which are not variable ac- 
cording to our investigations are marked by a cross (+). The 
two regions NGC 2264 N & S with the highest concentration 
of variables are marked by boxes. 



will be a certain degree of contamination by non-PMS stars 
if one selects candidates on the basis of variability. Therefore 
we have to disentangle variable PMS stars from variable back- 
ground or foreground stars. This was done in two steps, namely 
using first the Ic vs (Rc-Ic) colour-magnitude diagram (CMD) 
and second the (Rq-Ho) vs (Rc-Ic) colour-colour diagram. In 
the following subsections we describe the two discrimination 
procedures in detail. Variable stars which pass both of these 
tests will be likely PMS members of the cluster and for sake of 
simplicity will therefore be called PMS stars in the subsequent 
discussion. For a full confirmation of the PMS nature of these 
stars additional observations are necessary, like measurements 
of the Li I /16707-line equivalent width. We expect, however, 
that only a minor fraction of the stars which passed both of the 
tests are non-PMS stars. 



5.1. PMS test I: The colour-magnitude diagram and 
determination of the PMS region 

In a CMD the PMS stars are located in a region above the 
main sequence (MS). This is due to their larger stellar radii 
and correspondingly larger brightness compared to MS stars 
of the same spectral type or colour. In order to determine the 
borders of the PMS region we placed a well defined sample of 
PMS stars in the CMD. This sample was selected in such a way 
that it is contaminated as little as possible by background and 



foreground stars. It consists of two subsamples: 1) previously 
known cluster members selected from the catalogues of Park, 
Sung & Bessel d2000b and Sung, Bessel & Lee ( fT997b and 2) 
a well selected subsample of the newly found periodic and ir- 
regular variables. The second subsample consists of our newly 
found periodic and irregular variables located in the two dense 
concentrations of variables in the cluster. These two concen- 
trations are called NGC 2264 N & S and are marked in Fig.0 
The two regions are a part of the so called "on cloud region" 



defined by Rebull et al. ( 2002 1 and as already mentioned above 
are known as regions of maximum stellar density in the cluster 
(Mathieu [T986l Sagar et al. [T^8l . 

Determining the borders of the PMS region in the CMD 
only from the variable stars in NGC 2264 N & S has the advan- 
tage of producing a maximum fraction of PMS stars relative to 
any background stars. This is due to high extinction of the dust 
located behind the cluster stars. This is shown by Fig. II llwhich 
also illustrates that in particular for NGC 2264 S a very large 
background extinction is present. Since we selected only vari- 
able stars located in these regions the contamination of the this 
subsample with (non-variable) MS foreground stars is also neg- 
ligible. Therefore the selected variable stars in the two regions 
NGC 2264 S and N are most likely PMS cluster members. 

The CMD of all variable stars in NGC 2264 N & S and all 
previously known PMS stars in our whole NGC 2264 field is 
shown in Fig. [8] The spectral types marked at the top of the 
panel were determined by using photometric measurements of 
ZAMS stars with known spectral type. For stars with spec- 
tral types earlier than MO mesurements reported by Johnson 
( 1966) were used. For spectral types later than MO measure- 



ments by Leggett ( 1992) and Kirkpatrick & McCarthy ( 1994) 
were used. If necessary, transformation into the Cousins sys- 
tem were performed by applying the transformation equations 
by Bessel ( 1983). For the spectral types zero reddening is as- 
sumed. From this diagram we defined the PMS region in the 
CMD. The lower and upper borders of this region are indicated 
by the dashed lines. 

The defined PMS region will obviously eliminate MS stars 
which are at the same distance as NGC 2264 or further away. 
Such stars are located in the CMD at or below the ZAMS indi- 
cated in Fig.|8]and are therefore well outside the PMS region. 
Because of their larger apparent Ic magnitude foreground MS 
stars with a distance between ~ 300 pc and ~ 600 pc could 
be located in the PMS region. However, as already mentioned 
above the contamination with these stars is probable negligible. 

On the other hand it is not possible to discriminate between 
variable PMS stars and variable background giants using the 
PMS region in the CMD if the latter are at (large) distances 
which results in apparent magnitudes that shift them into the 
PMS region. K or M giants are typically between 3 and 10 
magnitudes brighter than PMS stars at the same distance and 
the same spectral type. Therefore background giants with a cor- 
responding smaller apparent Ic magnitude could be located in 
the PMS region of Fig. [8] Since the reddening vector is nearly 
parallel to the borders of the PMS region even a highly red- 
dened giant will stay in the PMS region. In addition several 
types of giants are variable (e. g. RRLyrae stars) and therefore 
it is necessary to eliminate these stars from our sample of vari- 
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Fig. 8. Determination of the PMS region in the CMD. Only the 
variables in NGC 2264 S and N (see Fig.[7} and all previously 
known PMS stars in our whole NGC 2264 field are plotted. The 
symbols are the same as in Fig. [7] The two dashed lines indicate 
the lower and upper borders of the PMS region. 

Fig. 9. The {Rq-Ho) vs (Rq-Ic) colour-colour diagram for all 
stars which passed the first PMS test. Circles and triangles rep- 
resent periodic and irregular variables respectively. Filled sym- 
bols indicate the stars which passed both PMS tests. Variables 
which were previously known as PMS are indicated by a ver- 
tical line. Non-variable previously known PMS stars are indi- 
cated by a cross. The solid line represents the fit of the locus 
of PMS and MS stars. Variable stars below the dashed line 
(open symbols) were rejected from the analysis because they 
are probable background giants. 



able stars using a different selection criterion. This will be done 
in Sect-IP 

Of the 543 periodic variables we found in our whole 
NGC 2264 field, 451 stars were located in this PMS region 
and only these stars were used for the further analysis. These 
numbers together with the few (only 3) periodic variables 
outside the PMS region in Fig. [8] confirm that the regions 
NGC 2264 S & N are very little contaminated with foreground 
or background MS stars compared to the rest of our field. That 
furthermore indicates that the above described procedure to de- 
fine the PMS region in the CMD provides relatively reliable 
borders. 

5.2. PMS test II: The (R c - Ha) vs (R c - I c ) 
colour-colour diagram 

As outlined in the previous section it is essential to elimi- 
nate background giants from the sample of variables we found. 



Fig. 10. Locations of all PMS variables from our survey which 
passed both tests including non-variable previously known 
PMS members. The symbols are the same as in Fig. [5] Also 
shown (indicated by the two boxes) are the two regions 
NGC 2264 N & S which we used to determine the PMS region 
in the CMD. The star in the northern box indicates the location 
of the bright star S Mon (V=4.7 mag). 

Fig. 11. Locations of all stars in the field that failed at least one 
of our PMS membership criteria (see Sect. |5J|. We also show 
the two regions which we used for determination of the PMS 
region in the CMD. 

Fig. 12. The Iq vs (Rq - Ic) colour-magnitude diagram for dif- 
ferent samples of stars. In (a) we show all 543 periodic vari- 
ables we found. The 405 periodic variables that passed our two 
PMS criteria are marked with filled circles while those peri- 
odic variables which did not pass at least one test are marked 
by open circles. The PMS zone is indicated by the dashed lines. 
The solid line represents the ZAMS sequence. Panel (b) shows 
the CMD for the all 484 variable stars. The 184 irregular vari- 
ables which fulfilled both of our PMS criteria are indicated by 
filled triangles. The 300 variable stars which failed at least one 
test are shown by open triangles. In panel (c) we show the CMD 
of the previously known PMS stars in our whole NGC 2264 
field. 

Therefore we used the location of the stars in the (Rc - Ha) vs 
(Rc - Ic) colour-colour diagram (Fig.|9|l as a second additional 
selection criterion for classifying a star as a PMS star; i. e. all 
stars which passed the first test have also to pass the second 
one to be classified as a PMS star. Many PMS stars, in partic- 
ular CTTSs, show large Ha emission. Thus their (Rc - Ha) 
colour is larger than the (Rc - Ha) colour of a MS star of the 
same spectral type. There also exists WTTSs or "naked" TTSs 
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with weak or no detectable Ha emission (e.g. see Appenzeller 
& Mundt, 1989). The locus of these latter stars is similar to the 
MS in the (Rc - Ha) vs (Rc - Ic) colour-colour diagram. 

On the other hand giants with spectral types later than ~K3 
have smaller (Rc - Ha) colours than MS stars of the same 
(Rc - Ic) colour. The differences in (Rq - Ha) between gi- 
ants and PMS/MS stars result from stronger molecular bands 
in the giants' spectra causing different spectral energy distribu- 
tions in the Rc band. In addition giants are typically at larger 
distances and are therefore highly reddened. Since reddening 
affects mainly the {Rq - Ic) colour they are shifted further away 
from the MS (i. e. to the right) with increasing distance. This 
behaviour was confirmed by a simulated (Rc-Ha) vs (Rq -Ic) 
colour-colour diagram with the help of standard star spectra of 
different spectral type and luminosity (i. e. by multiplying the 
filter transmission curves of the WFI with the spectral energy 
distribution of standard stars spectra). 

Therefore it is possible to discriminate between (back- 
ground) giants and PMS stars using their location in the (Rc - 
Ha) vs (Rc - Ic) colour-colour diagram. We determined the 
median (Rc - Ha) colour as a function of (Rc - Ic) for the lo- 
cus of PMS/MS stars using the data shown in Fig. [5] In Table0 
we list the median (Rc - Ha) colour as a function of (Rc - Ic) 
for the stars on the PMS/MS locus. These values were derived 
in different (Rc - Ic) colour bins of 0.05 magnitude width. We 
have fitted a quadratic function to these medians and obtained 
a locus relation of the form 

(Rc-Ha) locm = -0.06 x (R c -I c ) 2 + 0.38 x(R c -I c )- 3.50.(6) 

This relation is shown in Fig. [9] as the solid line. As we will 
show in Sect. 16. 21 . the vertical distance of a star from this line 
measures the strength of the star's Ha emission. 

In order to discriminate between PMS stars and background 
giants we defined a lower discrimination level shown in Fig. [9] 
as a dashed line. This line is given by (Rc - Ha)i ov/ - 1.65 x 6, 
where 5 is a fit of the standard deviation in (Rc - Ha) in each 
(Rc-Ic) bin and is given by 6 = e 0.60(« c -/c)-3.42 > Assuming that 
the locus of stars in each bin are Gaussian distributed 95% of 
the locus stars are located above the lower discrimination level 
which therefore represents a lower envelope of the locus. Any 
variable star (periodic or irregular) which passed the first se- 
lection criterion described in Sect. 15. H and that is located above 
this level was kept for the final analysis while stars below this 
level were rejected. 

From the 451 periodic variables that passed the first se- 
lection criterion another 46 were rejected because of their 
(Rc - Ha) colour. We finally have a total of 405 periodic and 
184 irregular variables (with P vfln > 99.9%) which are all very 
likely PMS members of the NGC 2264 star forming region. 
Hence, we call these stars variable PMS members. The peri- 
odic variable PMS members are separately marked in Table [6] 

5.3. The spatial distribution of PMS members 

In Fig. 1101 and Fig.^2 we show the spatial distribution of all 
PMS stars and the positions of non-cluster members, respec- 
tively. The latter were selected from the complete list of all 
10554 stars in the field and represents all stars that failed at 



least one of our PMS membership criteria described in the 
previous sections. From the spatial distribution of the non- 
cluster members one can clearly identify the extent of the dust 
cloud located towards NGC 2264. When comparing Fig. [TO] 
with Fig. [^it is evident that most PMS stars are located in 
the region with the highest background extinction, i. e. close to 
the dense gas and dust out of which they have probably been 
formed. 

5.4. Colour-magnitude diagrams of the periodic and 
irregular variables 

In Fig.[^]we show the Ic vs (Rc - Ic) colour-magnitude di- 
agram of all periodic and irregular variables. It is evident that 
the periodic variable PMS members are much more concen- 
trated around a line parallel to the ZAMS than the irregular 
variable PMS members. The reason for this could be a smaller 
age range of the periodic variables and/or a higher variability 
of the irregular variables compared to the periodic variables. In 
addition a higher intrinsic extinction of the irregular variables 
due to circumstellar disks could lead to a larger scatter. In the 
following section the differences between the two subsamples 
will be discussed in more detail. 

In Fig. ^] we also show the Ic vs (Rc - Ic) colour- 
magnitude diagram of all (169) previously known PMS stars 
in NGC 2264 which we could identify in our sample. It is ev- 
ident that the new PMS variables we found extend to much 
fainter magnitudes than the previously known PMS stars, i. e. 
we found new PMS stars mainly in the low mass regime which 
probably reaches down to the substellar limit of approximately 
Ic = 18.5 mag in NGC 2264. 

6. The nature of the variability 

In this section we investigate some general aspects of the stars' 
variability and investigate which differences in physical prop- 
erties between periodic and irregular variables can be found. 

6.1. The degree of variability 

We use the standard deviation cr of the light curves as an es- 
timation of the degree of variability. The advantage of using 
cr (see Sect. I4.2> rather than the peak-to-peak (ptp) variation 
of the stars for the characterisation of the variability is that cr 
is more robust to outliers in the light curves (e. g. due to cos- 
mic rays). In order to account for those statistical outliers in 
the light curves we have equated the ptp variation for each star 
with the magnitude difference between the third highest and 
the third lowest data point. Both quantities cr and ptp variation 
are given in Table 0] for each star. We find that the ptp vari- 
ation is typically three times the standard deviation but it can 
vary between two and five times the standard deviation in some 
extreme cases. 

In Fig. ^] we show the cumulative distribution of the stan- 
dard deviation cr for all 405 periodic and 184 irregular PMS 
variables. It is evident from this plot that about half of the pe- 
riodic variables have cr < 0.02 mag while only 5.4% (10/184) 
irregular variables show this small level of variability. On the 
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Fig. 13. (a) The top panel shows the cumulative distribution of 
the standard deviation cr for all 405 periodic and 184 irregular 
variable PMS stars respectively, (b) The lower panel shows the 
histogram of cr for the two different samples of stars in (a). 

other hand only 2.7% of the periodic variables have a cr > 

0. 2 mag but 14.7% of the irregular variables exceed this level, 

1. e. the fraction of stars with a small cr is always higher for the 
periodic variables. 

This excess of periodic low amplitude stars is not surpris- 
ing because the advantage of the periodogram analysis is that 
it is sensitive even for low amplitude variations. This results 
from the fact that the power from photometric errors (white 
noise) is distributed over the whole frequency domain of the 
power spectrum. Therefore even low amplitude variations close 
to (and below) the photometric errors can be detected. The x 1 
test on the other hand is able to detect variability only if the 
variations in a light curve are inconsistent with the photometric 
errors. Our cutoff level of a probability P var j = 99.9% corre- 
sponds to deviations of more than 3cr. Therefore we expect to 
find more periodic variables than irregular variables among the 
low amplitude variables. 

In order to take any bias resulting from the different vari- 
ability analysis methods into account we selected only those 
264 periodic variables which are also variable according to the 
X 2 test and compared their distribution with that of the irregular 
variables. We found again that the two distributions are signif- 
icantly different, as confirmed by a Kolmogorov-Smirnov test 
(Press et al. 1992) which shows that the probability that the 



two distributions are equivalent is less than 5 x 10~ 15 . A further 
discussion on this follows in Sects. l6~2l and l9l 

A similar result was found in the ONC (Herbst et al. 2002 1 
and in other T Tauri star associations. Herbst et al. 1 1994 1 have 
shown that most irregular variables are CTTSs and that they 
show large variations in their light curves. Periodic variables 
on the other hand are mainly WTTSs and have peak-to-peak 
variations in Ic of less than 0.5 mag. This would correspond to 
standard deviations of cr = 0. 1 to 0.25 mag for stars in our sam- 
ple. Herbst et al. (|2002i furthermore concluded that the varia- 
tions in the stellar brightness of the periodic WTTSs are mainly 
caused by cool spots and the observed maximum peak-to-peak 
variation of 0.5 mag is interpreted as the maximum possible 
brightness change which cool spots on the surface of a K or 
M star could cause. The periodic variables with peak-to-peak 
variations larger than 0.5 mag are believed to be due to hot sur- 
face spots resulting from mass accretion. In our sample 3.5% of 
the periodic variables have peak-to-peak variations > 0.5 mag 
(cr > 0.16 mag). The irregular variability of CTTSs can have 
various reasons. Often it is attributed to variable mass accretion 
resulting in hot spots which are not stable in brightness, size, 
and location over a few rotation periods. In addition flare-like 
activities can be an additional source of variability. In reality 
the situation is probably more complicated and it could well be 
that some WTTSs have small hot spots and that CTTSs have 
cool spots in addition to large hot spots. However, for CTTSs 
it is much more difficult and in many cases even impossible 
to detect periodic brightness modulations caused by cold spots 
because of the overlying "noise" from irregular variability. 

6.2. The Ha emission index of periodic and irregular 
variables 

If the fraction of CTTSs is in fact higher among the irregular 
variables they should also show larger Ha emission than the 
periodic variables. In order to investigate the Ha emission of 
the stars in the two samples we defined a Ha emission index 
A(Rq - Ha) of a star by the following equation: 



A(R C - Ha) = (R c - Ha) stm - (R c - Ha)^, 



(7) 



where the (Rc - Ha)\ ocus is shown as a solid line in Fig.|9]and 
is given by Eq. |6] It is the (Rc - Ha), (Rc - Ic) relation of a 
MS star. The Ha index is a measurement of the Ha emission. 
Using stars with known Ha equivalent width (W A (Ha)) we find 
(see Paper II) that 83% of the stars with an emission index of 
A(R C - Ha) > 0.1 mag have W A (Ha) 10A and are therefore 
likely CTTSs. On the other hand about 85% of the stars with 
A(R C - Ha) < 0.1 mag have W A (Ha) < 10A and are therefore 
most likely WTTS. 

In Fig.^^we show the cumulative distribution and the his- 
togram of the Ha index for periodic and irregular variables. 
It is evident that the Ha emission index of the periodic vari- 
ables is concentrated around zero and that only 22% of these 
stars exceed an Ha index of 0.1 mag. On the other hand the 
Ha emission index of the irregular variables is distributed over 
a larger range and for 68% of the irregular variables it is above 
the critical level of 0.1 mag. 
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Fig. 14. (a) The top panel shows the cumulative frequency dis- 
tribution of the Ha index A(Rc - Ha) for all 405 periodic and 
184 irregular variable PMS stars, (b) The lower panel shows 
the histogram of the Ha index for the two different samples of 
stars in (a). 
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Fig. 15. The Ha index as a function of the degree of variabil- 
ity measured by the standard deviation cr for all 405 periodic 
and 184 irregular variable PMS stars. Circles represent periodic 
variable stars while triangles represents irregular variables. The 
dashed and solid lines represent the medians of the Ha index 
for the irregular and periodic variables respectively in <x-bins 
of variable width with each containing 14 data points. 



This supports the interpretation that the fraction of CTTSs 
among the irregular variables is higher than the fraction of 
CTTSs among the periodic variables and vice versa for the 
WTTSs. In Sect. [9] we will quantify how strongly the sample 
of periodic variables is biased towards the WTTSs by compar- 
ing the fractions of CTTSs and WTTSs among the periodic 
variables with the expected fractions. 

6.3. The correlation between cr and the Ha index 

In order to investigate this interpretation further we look for 
a correlation between the degree of variability of the stars 
and their Ha index which is a measurement of their accretion 
and/or chromospheric activity. In Fig.^l we show the Ha in- 
dex and its median as a function of cr for periodic and irregular 
variables respectively. It shows that for both samples there is 
an increase of the median Ha index with increasing cr. Each 
median was calculated in cr bins of variable width containing 
14 data points. Only the range cr < 0.22 mag is considered here 
since there are insufficient data points for larger values of cr. 
A Spearman rank-order correlation test (Press et al. 1992jl in- 
dicates that the probability that cr and the Ha index are not 
correlated is less than 6 x 10~ 28 if we use the joined sample of 



periodic and irregular variable stars. However, this high prob- 
ability that there is a correlation between the two quantities is 
dominated by the periodic sample. If we use only the periodic 
sample the probability that cr and the Ha index are not corre- 
lated is less than 9 x 10~ 15 . On the other hand the probability 
that cr and the Ha index are correlated for the irregular vari- 
ables is only 0.863. 

It is interesting that variability in the periodic variables is 
strongly correlated with the Ha index while there is no evi- 
dence for a correlation for the irregular variables. One possible 
reason is that the sizes and numbers of cool spots (i. e. the total 
spot area) on periodic variables (WTTSs) are correlated with 
the chromospheric activity as in other active late type stars, i. e. 
periodic variables with more spots have more active regions on 
their surface and therefore stronger chromospheric Ha emis- 
sion. For the irregular variables we see that there is a large 
scatter in the Ha, cr relation, and even low amplitude variables 
have large Ha indices. This could imply that there are many 
CTTSs which have no large variations in their mass accretion 
rates. In addition these stars must have hot spot patterns which 
do not cause any large light modulations (e. g. accretion rings 
symmetric to the rotation axis, see e. g. Mahdavi & Kenyon 
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Fig. 16. The (V - Rc) vs (Rc - Ic) colour-colour-diagram for 
periodic (top panel) and irregular variables (bottom panel) that 
passed both of our two PMS tests. Those stars which have an 
Ha-inde\> 0.1 and therefore show the properties of accreting 
CTTSs are marked with a surrounding circle. The solid lines in 
both panels represents the ZAMS. The arrow length and direc- 
tion indicates the mean reddening towards NGC 2264. 



6.4. The (V - R c ) vs (Rc - Ic) colour-colour-diagram 

In Fig. we show the (V - Rc) vs {Rc - Ic) colour-colour- 
diagram for periodic and irregular variables which passed both 
of our two PMS tests. The stars with an Ha-index of A(Rc - 
Ha) < 0.1 are marked separately. It is evident that the scat- 
ter is much larger among the irregular variables compared to 
the periodic variables. Agreement with the ZAMS is good, al- 
though deviations due to the smaller log g values of the PMS 
stars are expected. The data points below the MS (i. e. with 
larger (V - Rc) colours than a MS star) can be explained in 
both plots by embedded stars which are highly reddened with 
up to eight times the mean reddening. This interpretation is sup- 
ported by the analysis of Rebull et al. (2002 1. They have deter- 
mined an average E(R C - I c ) = 0.1 + 0.02 mag in NGC 2264, 
but showed that some stars in NGC 2264 are reddened by up to 
E(Rc - Ic) = 1 mag. 
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Fig. 17. The (Rc - Ha) vs (Rc - Ic) colour-colour-diagram 
for all "non-variable" stars which passed the first PMS test 
(see Sect. including previously known "non-variable" PMS 
members. The dotted line represents the modified second se- 
lection criterion for "non-variable" PMS stars (see text). All 
"non-variable" stars which are located by more than their error 
in (Rc - Ha) above this level are classified as PMS candidates. 



7. "Non-variable" PMS members with strong Ha 
emission 

Based on our periodic and irregular variability study we have 
selected a total number of 584 PMS members (405 periodic 
and 184 irregular variables) using the two selection criteria 
described in Sect. |5] However, there may exist low-amplitude 
PMS stars, for which we could detect no significant variability. 
Reasons for a small degree of variability in WTTSs could be a 
more or less equally distributed pattern of many spots causing 
only a small modulation of the stellar brightness and in CTTSs 
there could be phases of relatively stable mass accretion. To 
simplify matters we call these stars "non-variable" PMS stars. 
In order to search for "non-variable" PMS stars both PMS tests 
described in Sect.[5]were applied to all stars which were not de- 
tected to be variable. "Non-variable" PMS stars were selected 
if they passed both tests and have in addition strong Ha emis- 
sion. 

In Fig.ll7lwe show the (Rc-Ha) vs (Rc~Ic) colour-colour- 
diagram for all of our non-variable stars which passed the 
first selection criterion. From these stars we selected only stars 
which showed an enhanced Ha emission using a more conser- 
vative threshold in the colour-colour diagram. This threshold 
is indicated by the dotted line in Fig.[n]which represents the 
upper envelope of the MS locus (solid line). The dotted line is 
given by 



(R c - Ha) U gh = (Rc - tfaOiocus + 1.65 x 6, 



(8) 



M. H. Lamm et al.: A Rotational & Variability Study of PMS Stars in NGC 2264 



(for the definition of 5 and (Rc - Ha)\ ocus see Sect. 15. 2> . This 
upper discrimination level is equal to the lower discrimination 
level used in Sect. l5.2l mirrored at the PMS locus. 

We consider only stars as "non-variable" PMS members 
if they are located more than their photometric error A(R C -Ha) 
above the upper discrimination line, i. e. (Rc- Ha)- A(R C -Ha) ^ 
(Rc - Ha)hi g h- In total 32 stars of the non-variables with 
(Rc - Ic) < 2.5 mag in Fig.[T7]passed this additional PMS test. 
Out of these "non-variable" PMS members 2 stars (no. 5941 
and no. 6164) were classified as PMS stars prior to this study. 
Another 12 previously known "non-variable" PMS members 
did not pass this additional test. The 30 newly found stars are 
separately marked in Fig.ll7land listed in Table [S] 

As one can see from Fig.[n]most (83%) of the new "non- 
variable" PMS members have (Rc -Ic) ^ 1.8 mag. This can be 
understood if we consider that the photometric errors for these 
stars are typically larger compared to bluer (i. e. brighter) stars. 
Therefore any variations in the light curve are harder to detect. 

In total we have identified 621 PMS stars in NGC 2264 
including 405 periodic variable 184 irregular variable and 
32 "non-variable" stars with relatively strong Ha emission. 
Compared to the 182 previously known PMS stars and PMS 
candidates reported by Park et al. 2000 (not including 26 previ- 
ously known massive OB stars) we have increased the number 
ofknownPMS stars in NGC 2264 by a factor of 3.4. As already 
mentioned in Sect. !5.3l most of the newly found PMS stars are 
low mass stars (M 55 0.25 Mq) with masses probably extending 
into the substellar regime. 

8. Fraction of variable PMS stars and 
completeness level of our PMS sample 

In this section we estimate the fraction of variable PMS stars 
in the cluster. From this fraction we estimate the complete- 
ness level of our sample of variable PMS stars in NGC 2264. 
Therefore the central question in this section is: what fraction 
of PMS stars did we find through our photometric monitoring 
program and the adopted methods in selecting PMS stars? 

The fraction of variable cluster members is given by / var ; = 
Nvari/Wvari + ^V non _ V aii), where Af vai i is the number of variable 
(both periodic and irregular) and A^ n0 n-vari is the number of 
"non-variable" PMS stars in NGC 2264. These numbers (in 
particular A^on-van) are well known for CTTSs since the PMS 
nature of these stars can be confirmed by an enhanced Ha emis- 
sion. In addition we expect that there are also WTTSs (with no 
or weak Ha emission) which are "non-variable"; i. e. they have 
light modulations below our detection limit. Using the Ic vs 
(Rc - Ic) colour-magnitude or (Rc - Ha) vs (Rc - Ic) colour- 
colour diagrams these stars are indistinguishable from MS fore- 
ground stars. In order to estimate the number of "non-variable" 
WTTSs we therefore have to keep the relative contamination 
by foreground MS stars as low as possible. This is achieved by 
using only the two regions NGC 2264 N and S (see Sect. 15. 2> 
for our analysis rather than the whole observed field. In or- 
der to be not affected in the following analysis by photometric 
errors we restrict ourselves to stars with Ic < 18.0 mag and 
(Rc-Ic) < 1.8 mag. 
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Table 9. The estimated fraction of variable stars in 
NGC 2264 N and S with I c < 18.0 mag and (R c - I c ) < 
2.0 mag. Listed are number and fractions of periodic PMS vari- 
ables (A^p_vari), irregular PMS variables (M-varO, "non-variable" 
PMS candidates (A^on-vanX the total number of all PMS candi- 
dates (NpMS-totai = Np-vari + M-vari + Mion-vari), and the number 
and fractions of PMS variables (Np-vari +M-van) in the two con- 
centrations of PMS stars NGC 2264 N and S (for a definition of 
the PMS candidates see text). The fractions were calculated rel- 
ative to the total number of PMS star candidates (A^PMS-totai) in 
the two regions. 



sample of stars 


North 


South 


North & South 


-NpMS-total .... 
N 

JY non-van 

/V 

M-vari 


110 


120 


230 


25 (22.7%) 
56 (50.9%) 
29 (26.4%) 


36 (30.0%) 
54 (45.0%) 
30 (25.0%) 


61 (26.5%) 
110(47.8%) 
59 (25.7%) 


variables .... 


85 (77.3%) 


84 (70.0%) 


191 (73.5%) 



In Table[9]we list the number of variable and "non-variable" 
stars in each of the two regions which passed both PMS tests 
(see Sect. and are in these colour and magnitude ranges. 
From the last column of this table it is evident that in the two 
regions NGC 2264 N and S about 74% of the PMS candidates 
(i. e. stars that passed the two tests) are variable. About half 
of the PMS candidates show periodic light modulations while 
only one fourth of the PMS candidates are irregular variables. 
We note that these fractions are only lower limits since we as- 
sume that all stars in the two regions which passed both tests 
are indeed PMS stars. In particular we assume that the contam- 
ination by non-members in the variable sample is negligible; 
i. e. all variables that passed both PMS test are actually PMS 
stars. 

However, if we further assume that the two regions are rep- 
resentative for the whole cluster we can conclude that at least 
/vari = 74% of the PMS stars have been detected to be variable. 
This lower limit of the fraction of detected variable stars is not 
significantly higher if we consider only brighter stars; e. g. for 
stars with I c < 16.0 mag in NGC 2264 N and S we obtain a 
fraction of / var i =76%. 

Since we expect all PMS stars in the cluster to be variable at 
some level our variability study is complete at a 74% level for 
stars with Ic < 18.0 mag and (Rc - Ic) < 1.8 mag. As we will 
show in the following section the completeness level differs for 
WTTSs and CTTSs. While nearly all (95%) of the CTTSs are 
found to be variable only 68% of the WTTSs are found to be 
variable. 

In summary we conclude that our method of photometric 
monitoring is a very powerful tool for finding most PMS stars 
in our cluster since at least three quarters of the whole PMS 
population in NGC 2264 (with (R c - Ic) < 1.8 mag) could be 
identified this way. Therefore, our database of periodic and ir- 
regular variables should be a representative subset of the PMS 
stars in NGC 2264. 
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9. Fractions of periodic variables among CTTSs 
and WTTSs 

In the following we try to estimate the fraction of periodic vari- 
ables among the WTTSs and CTTSs in order to investigate 
how strongly the period distribution of NGC 2264 is biased 
by the non-detection of periods in these two groups of stars. 
Again we restrict our analysis to stars with Iq < 18.0 mag and 
(Rc - Ic) < 1-8 mag. Furthermore we consider in this section 
only stars that passed both PMS tests. 

As already outlined in Sect.[6]the periodic variables are bi- 
ased towards the WTTSs since it is much harder to detect peri- 
odic brightness modulations of CTTSs in the presence of a su- 
perimposed irregular variability. However, we will show below 
that this bias towards the WTTSs will be partially compensated 
by those WTTSs which have brightness modulations below our 
detection limit and have therefore not been found. 

As in Sects. |6]and0we regard all stars with enhanced Ha 
emission measured by the Ha index as CTTSs. This is the 
case for A^ctts = 145 stars in our whole investigated region. 
While 89 (61.4%) of these stars are irregularly variable only 
49 (33.8%) are periodically variable and 7 (4.8%) are "non- 
variable". Thus only for one third (49/145) of the CTTSs peri- 
ods are detectable. 

It is more difficult to estimate the total number of WTTSs 
in the cluster (AVtts) since "non-variables" WTTSs are indis- 
tinguishable from MS foreground stars in the Ic vs (Rc - Ic) 
colour-magnitude and in the (Rc - Ha) vs (Rc - Ic) colour- 
colour diagram. However, with the results of the previous sec- 
tion we can relatively easily estimate the total number of PMS 
stars in the cluster (Npms) and since A^pms = Afcrrs + AVtts we 
can determine A^wtts indirectly. 

Therefore we first estimate A^pms- In total 465 stars with 
(Rc-Ic) ^ 1.8 mag and /c S 18.0 mag are periodic or irregular 
variables. Since at least 73.5% of the PMS stars are variable 
(see Sect. Hi we get N PM s = 465/0.735 - 630. With AVtts = 
A^pms - A^ctts we g et AVtts = 630 - 145 = 485 which is an 
upper limit since 0.735 is a lower limit. 

Of the periodic variables with (Rc - Ic) ^ 1.8 mag and 
Ic < 18.0 mag 265 have a small Ha index and are therefore 
classified as WTTSs. In addition 62 irregular variables are clas- 
sified as WTTSs. Using these numbers we conclude that for at 
least 55% (265/485) of the WTTSs in NGC 2264 we have de- 
tected a period and at least 13% (62/485) of the WTTSs are 
irregular variable; i. e. for 68% of the WTTSs we detected vari- 
ability. 

How do the measured period fractions of 55% and 34% 
for the WTTSs and CTTSs respectively influence the final pe- 
riod distribution? To answer this question let us first assume 
that the periods for all (estimated) 630 PMS stars in the cluster 
with (Rc - Ic) ^ 1.8 mag and Ic < 18.0mag could be mea- 
sured. Using Ncms = 145 and AVtts = 485 (see above) we 
get that 23% (145/630) of the stars are CTTSs while the frac- 
tion of WTTSs is 77% (485/630). In our study 16% (49/341) 
of the periodic variables in these colour and magnitude ranges 
are CTTSs and 84% (265/314) of the periodic variables are 
WTTSs. From these numbers we conclude that the final period 
distribution is slightly biased towards WTTSs. 



10. Summary & Conclusions 

We have carried out an extensive search for rotation periods and 
variability in NGC 2264. The main results of this investigation 
are as follows: 

1. We have obtained absolute photometry in the Cousins 
V, Rc, and Ic bands and instrumental Ha magnitudes 
for about 10600 stars with 9.8 mag < Ic < 21 mag in 
NGC 2264. Relative light curves with more than 80 data 
points each were obtained for each of these stars in the Ic 
band. Two different periodogram analysis techniques and a 
X 1 test yielded a sample of 543 periodic and 484 irregular 
variables with 11.4 mag < Ic < 19.7 mag in the observed 
field. 

2. In order to check the PMS nature of the periodic and ir- 
regular variables two different PMS selection criteria were 
applied. These criteria are their locations in the Ic vs 
(Rc~Ic) colour-magnitude diagram and their location in the 
(Rc - Ha) vs (Rc _ Ic) colour-colour diagram. In this way 
405 periodic and 184 irregular variable PMS stars could be 
selected from the sample. This is a enormous increase in the 
number of known rotation periods compared to the about 30 
published periods prior to this study. The rotation periods 
of the periodic variables are typically between a quarter of 
a day and 15 days. 

3. Using an enhanced Ha emission index as a PMS indicator 
we could identify an additional 35 PMS stars which are 
"non-variable" according to our study. Together with the 
589 variable PMS stars we found a total of 624 PMS stars 
in NGC 2264. Before our study only 182 PMS members 
and candidates were known. Most of our newly found PMS 
stars are fainter than Ic - 15 mag and of late spectral type 
(£M2). 

4. Comparison of the distribution of the standard deviation cr 
in the light curves as well as the distribution of the Ha 
emission index for periodic variables with the correspond- 
ing distributions for the irregular variables strongly sug- 
gests that our method for determining the rotation peri- 
ods preferentially selects weak-line T Tauri stars (WTTSs) 
while classical T Tauri stars (CTTSs) are underrepresented. 
This interpretation is confirmed by the calculated fractions 
of periodic variables among WTTSs and CTTSs which, 
however, show that the bias of the measured rotation peri- 
ods towards WTTSs is relatively small. We have estimated 
that 23% of the PMS stars in the cluster are CTTSs but only 
15% of the periodic variables are CTTSs. The correspond- 
ing numbers for WTTS are 77% and 85% respectively. 

5. The different properties of the two subclasses of T Tauri 
stars (CTTSs and WTTSs) are supposed to be due to the dif- 
ferent mechanisms which produce the variability. While the 
variability of WTTSs is likely due to large (magnetically) 
cool spots the variability of CTTSs is believed to be mainly 
caused by hot spots resulting from accretion flows onto the 
star. On CTTSs the periodic brightness variations resulting 
from cool spots (which are believed to exist on their surface 
too) are often undetectable due to the "noise" from the ir- 
regular variability. Determining the rotational periods for a 
larger fraction of CTTSs by photometric monitoring would 
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therefore require more extensive observations (e. g. larger 
time base and dense sampling). 
6. We have estimated that we have detected variability among 
at least 70% of the PMS stars (defined according to our 
PMS tests) in NGC 2264 with I c < 18.0 mag and (R c - 
Ic) < 1.8 mag. From these numbers it is evident that exten- 
sive photometric monitoring with a similar accuracy as in 
our study is a powerful tool for finding most PMS stars. We 
emphasise that this method is highly efficient only if there 
is enough data with a proper sampling to allow a search for 
periodic variables because many PMS members were iden- 
tified by their periodicity and not just by their variability. 
We conclude that the sample of variable PMS stars obtained 
in this study is a representative subset of the cluster mem- 
bers and it is not expected that additional monitoring pro- 
grams will substantial increase the number of known PMS 
stars at least for stars with Ic < 18.0 mag. 
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Appendix A: The CLEAN algorithm 

The CLEAN algorithm we used in Section 1 .21 subtracts the 
window function iteratively from the observed, or raw, spec- 
trum in the following way. First the spectral window function 



W(v) 



1 N 
N 



—2mvt„ 



and the observed spectrum 



D(v) — — d(t n )e 



-2mV/„ 



(A.l) 



(A.2) 



n=l 



are calculated at discrete frequencies with a spacing of Av = 
1/(10 T), where T = fi - is the total length of the data 
span. In our case T is about 62.8 days and the resolution of 
our CLEAN periodograms is therefore Av = 0.0016/<i. 

The input into the jth iteration is the spectrum Z),, where 
the iteration process starts with the (complex) raw spectrum 
Di = i = D. During each iteration step i = 1,2, ...,N the al- 
gorithm determines the frequency v, and the power |D;(v;)| of 
the highest peak in the spectrum £),-. Assuming that this peak 
represents a true period in the signal the aliasing powers A, at 
other frequencies are calculated with 



A,(v) = Ci W(v - vd + (ciTW(y + v ; ) 



(A.3) 



where c, is the /th complex clean component, that represents 
the amplitude of the peak and is defined as 



Di(yd - D;(v,)W(2v,) 



(A.4) 



l-|W(2v,)| 2 

(see Eq. (23) and (25) in Roberts et al. 1987 1. The c, were used 
in Eq. lA.3l for the calculation of the aliasing powers rather than 
the power |A(v,)| of the peak since the peak itself is affected 



by aliasing of the corresponding negative peak located at — v,-. 
The A, are subtracted from the spectrum D, and the resulting 
spectrum D,+i = D,- - A, is the input in the next iteration step. 

The factor g in Eq. (IA.4> is called the gain factor and takes 
care of the possibility that the highest peak itself could be an ar- 
tifact. Therefore only a fraction (g) of the power is used to cal- 
culate the alias powers. The iteration process stops after the Mh 
iteration if the power of the highest peak in the subtracted spec- 
trum Djv+i is below a given level P ml „ or if a maximum num- 
ber Nmax of iteration steps is exceeded. The last output Dn + \ 
is called the residual spectrum. We used the values g = 0.01, 
Pmin — 0.001 and N max = 300. Normally the peak power in the 
residual spectrum falls below this limit of P m ,„ = 0.001 before 
the maximal iteration of N max = 300 is reached. 

During the iteration process any artifact should be removed 
and the main output is the clean component spectrum which is 
defined as 



C(v) 



N r+ 

■ , xj — D> 



Ci 6(v - Vj) dv, 



(A.5) 



where 6(x) is the Dirac delta function and the c, are the clean 
components from Eq. (IA.4> . The resolution of the clean com- 
ponent spectrum in Eq. (IA.5> is given by the artificial frequency 
spacing Av used in the algorithm (see above). Since the reso- 
lution of the spectrum F that results from a discrete data set 
is Av ^ 1 /T we have to smooth the clean component spectrum 
with a beam function B(v). This is necessary because the power 
of the peaks in the spectra D[ is distributed into «b = Av/T fre- 
quency bins {points per beam) in the interval [v,— 1/T, v,+l/T]. 
The function B(v) is found by a fit of the main peak of the win- 
dow function at v = 0. 

The final step of the CLEAN algorithm is the calculation 
of the clean spectrum S(v) that contains the periodogram of 
the signal where the window function W(v) has been removed. 
S(v) represents the spectrum F(v) in Eq. Q. S(v) is obtained 
by adding the residual spectrum Dn+i(v) and the smoothed 
clean component spectrum C(v): 



S(v) = ZWv) + C(v). 
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Table 6. All periodic variables we found in our field. Listed are the identification number, the mean Ic magnitude, and the 
different periods detected by the Scargle periodogram technique (Ps e ), the CLEAN periodogram technique (Pcl) and by Kearns 
& Herbst ( 1998). Also shown is the period which we adopted in the end and the corresponding error for this period. Those 
periodic variables that passed both of our PMS criteria are marked with a "y" in the column "PMS member". Those stars that 
failed at least one of the two tests are marked with a "n" in the appropriate column. In the latter case the superscript indicates 
which test the star did not pass: a) rejected because of the location in Ic vs (Rc - Ic) diagram and b) rejected because of the 
(Rc - Ha) colour. The full table is available electronically. 



Star 


Ic 


(err) 




P CL P KH 


adopted P 


6P 


PMS 


6001 


16.91 


0.01 


6.01 


5.98 


6.01 


0.50 


n" 


6008 


16.70 


0.01 


1.32 


1.32 


1.32 


0.02 


y 


6012 


18.02 


0.03 


4.35 


4.33 


4.35 


0.26 


n'< 


6013 


14.57 


0.02 


10.48 


10.46 


10.48 


1.43 


y 


6014 


17.20 


0.01 


1.59 


1.59 


1.59 


0.03 


y 


6019 


16.67 


0.01 


0.65 


0.64 


0.65 


0.01 


y 


6022 


12.91 


0.02 


1.73 


1.73 


1.73 


0.05 


y 


6024 


15.38 


0.03 


9.71 


9.81 


9.71 


1.23 


y 


6031 


13.85 


0.01 


3.92 


3.90 


3.92 


0.20 


y 


6032 


14.72 


0.01 


0.80 


0.80 


0.80 


0.01 


y 


6039 


14.44 


0.01 


5.83 


5.76 


5.83 


0.44 


y 


6042 


16.91 


0.01 


0.54 


0.54 


0.54 


0.00 


y 


6043 


15.01 


0.01 


10.77 


10.64 


10.77 


1.51 


y 


6045 


17.19 


0.01 


1.75 


1.76 


1.75 


0.04 


y 


6055 


13.64 


0.01 


5.92 


5.98 


5.92 


0.47 


y 


6063 


15.33 


0.01 


3.38 


3.38 


3.38 


0.15 


y 


6064 


16.13 


0.01 


2.18 


2.17 


2.18 


0.07 


y 


6067 


12.78 


0.01 


2.93 


2.92 


2.93 


0.14 


y 


6077 


19.77 


0.03 


2.40 


2.37 


2.40 


0.08 


y 


6079 


15.83 


0.01 


5.22 


5.23 


5.22 


0.35 


y 


6081 


15.07 


0.01 


0.68 


0.68 


0.68 


0.01 


y 


6093 


16.53 


0.01 


0.96 


0.96 


0.96 


0.01 


n a,b 


6101 


18.00 


0.01 


0.46 


0.46 


0.46 


0.00 


n" 


6102 


15.15 


0.01 


9.04 


9.10 


9.04 


1.06 


y 


6115 


17.08 


0.01 


0.71 


0.71 


0.71 


0.01 


y 



Table 7. The relation between (Rc - Ic) and (Rc — Ha) colours for PMS/MS locus stars. The values for (Rc - Ha) are calculated 
in bins of 0.05 magnitude width, cr is the standard deviation of the (Rc - Ha) colours in each bin. 



(R C - Ic ) 


(R c - Ha) 


cr 


(R C - Ic ) 


(Rc ~ Ha) 


cr 


(Rc - Ic ) 


(Rc - Ha) 


cr 


(Rc - Ic) 


(Rc - Ha) 


cr 


0.38 


-3.36 


0.05 


0.83 


-3.22 


0.06 


1.28 


-3.10 


0.08 


1.73 


-3.02 


0.12 


0.43 


-3.33 


0.04 


0.88 


-3.24 


0.05 


1.33 


-3.10 


0.07 


1.78 


-3.03 


0.13 


0.48 


-3.33 


0.05 


0.93 


-3.22 


0.06 


1.38 


-3.08 


0.08 


1.83 


-2.99 


0.15 


0.53 


-3.34 


0.05 


0.98 


-3.22 


0.05 


1.43 


-3.09 


0.07 


1.88 


-3.02 


0.11 


0.58 


-3.30 


0.04 


1.03 


-3.20 


0.06 


1.48 


-3.07 


0.10 


1.93 


-3.02 


0.14 


0.63 


-3.30 


0.06 


1.08 


-3.17 


0.06 


1.53 


-3.07 


0.10 


1.98 


-2.92 


0.15 


0.68 


-3.30 


0.05 


1.13 


-3.16 


0.07 


1.58 


-3.06 


0.10 


2.03 


-3.03 


0.11 


0.73 


-3.26 


0.04 


1.18 


-3.13 


0.07 


1.63 


-3.04 


0.10 


2.08 


-2.89 


0.13 


0.78 


-3.24 


0.04 


1.23 


-3.12 


0.08 


1.68 


-3.05 


0.10 


2.13 


-2.87 


0.20 
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Table 8. New "non-variable" PMS stars in NGC 2264 selected using the PMS-test I and a slightly modified PMS-test II (see text). P vari is the p 
variable according to our^ 2 test. The other columns are the same as in Table|4] 



Star 


a( 72000") 

vj.\%j v j vj vj J 


6( 72000) 




err 


V - Ir 


err 


R 1 

i\q Ir 


crr 


iv L nil 


err 






p 

1 van 


cn 


3130 


6-40-21 32 

VJ a ^ VJ . Z. ± . _ 


9:48:26.5 


19.00 


0.03 






2.26 


0.03 


-2.31 


0.10 


0.040 


0.17 


0.0015 




3455 


6-40-26 95 

vj.^vj. j—\j.y^> 


9:49:55.7 


18.28 


0.02 






2.19 


0.02 


-2.57 


0.05 


0.014 


0.05 


32 81 53 

Ji. U1JJ 




3646 


6:40:29.41 


9:47:37.5 


16.83 


0.02 


3.26 


0.06 


2.07 


0.02 


-2.13 


0.07 


0.017 


0.05 


99.3415 


r- 


3648 


6-40-29 42 

VJ . I VJ . V . 1 


9:41:04.4 


18.60 


0.01 


3.39 


0.07 


2.07 


0.03 


-2.53 


0.06 


0.017 


0.06 


0005 

vj .vjvjvj^j 




3680 


6-40-29 89 


9-50-20 2 


16.79 


0.01 


2.62 


0.01 


1.50 


0.01 


-2.90 


0.01 


0.007 


0.03 


0000 

VJ ■ VJVJVJ VJ 




3840 


6-40-31 92 

vj • ™vj • Ji.yz. 


9:50:16.2 


19.42 


0.03 






2.28 


0.05 


-2.54 


0.15 


0.034 


0.12 


1.6468 




4323 


6-40-38 30 

VJ • ™VJ • JO. JU 


9:29:25.8 


17.71 


0.04 


3.06 


0.02 


1.81 


0.01 


-2.76 


0.05 


0.157 


0.65 


0363 

VJ . VJ J VJ J 




4877 


6-40-45 89 

vj « rvj > uVj y 


9:39:19.6 


19.19 


0.04 


2.07 


0.08 


2.09 


0.09 


-2.53 


0.18 


0.035 


0.13 


0.0031 




4924 


6-40-46 68 

VJ . "T VJ . T^VJ . VJ (J 


9:48:13.0 


16.92 


0.01 


3.05 


0.02 


1.86 


0.02 


-2.46 


0.02 


0.007 


0.03 


0000 

VJ .VJVJVJVJ 


Rl 


4960 


6:40:47.17 


9:28:50.4 


11.45 


0.01 


0.88 


0.01 


0.44 


0.01 


-3.10 


0.06 


0.014 


0.06 


0000 

VJ • VJVJVJVJ 


y; 


5456 


6-40-54 74 

VJ . ^ V J a -J ^ ■ / 


9:54:57.7 


19.23 


0.02 






2.23 


0.07 


-1.95 


0.44 


0.034 


0.17 


0.7261 




5567 


6-40-56 64 

VJ ■ ™VJ • *J VJ . VJT^ 


9-32-20 2 


19.33 


0.03 






2.30 


0.05 


-2.34 


0.16 


0.024 


0.10 


0000 

VJ . VJVJVJVJ 




5570 


6-40-56 73 

VJ . ^VJ . .J VJ . 1 .J 


9:38:10.7 


18.98 


0.02 






2.48 


0.07 


-1.96 


0.10 


0.031 


0.14 


59.4923 




5627 


6-40-59 39 

VJ ■ t^vj .^J J .J J 


9:31:00.2 


17.60 


0.02 


3.14 


0.09 


1.86 


0.02 


-2.32 


0.03 


0.022 


0.07 


99.3212 




5631 


6-40-59 53 

VJ . -tvj . j J . U *J 


9:34:44.2 


19.23 


0.02 






2.22 


0.05 


-2.32 


0.09 


0.044 


0.16 


3.4192 




5982 


6-41-07 56 

VJ . ^ L . VJ / .^JVJ 


9:41:34.1 


16.07 


0.01 


2.53 


0.04 


1.46 


0.01 


-2.51 


0.10 


0.014 


0.05 


2.5542 


R: 


6047 


6-41-09 23 

VJ . ^ ± . VJ V ■ .J 


9:29:05.0 


18.48 


0.02 






2.11 


0.03 


-1.99 


0.06 


0.019 


0.07 


0.2423 




6169 


6:41:13.11 


9:24:36.9 


16.44 


0.01 


2.79 


0.01 


1.61 


0.01 


-2.88 


0.01 


0.006 


0.03 


0000 

VJ .VJVJVJVJ 


Y_" 


6190 


6:41:13.77 


9:29:32.5 


17.64 


0.02 


3.53 


0.04 


2.12 


0.02 


-2.65 


0.03 


0.009 


0.04 


0.0394 


R: 


6306 


6-41-17 38 


9-35-56 5 


19.28 


0.03 






2.44 


0.05 


-2.52 


0.12 


0.030 


0.13 


0.0551 




6354 


6-41 -18 66 

VJ.T^ ± . A U . VJVJ 


9:49:02.7 


18.20 


0.01 


3.31 


0.04 


1.92 


0.01 


-2.73 


0.06 


0.017 


0.06 


0.4781 




6355 


6:41:18.69 


9:32:52.7 


19.03 


0.02 






2.32 


0.07 


-2.25 


0.08 


0.026 


0.12 


0000 

VJ .VJVJVJVJ 




6359 


6:41:18.77 


9:35:19.6 


19.45 


0.04 






2.33 


0.07 


-2.47 


0.16 


0.033 


0.14 


0.3219 




6490 


6:41:21.44 


9:53:01.9 


17.86 


0.01 


2.94 


0.03 


1.72 


0.03 


-1.81 


0.04 


0.018 


0.07 


43.0053 




6607 


6:41:23.29 


9:32:30.7 


18.70 


0.03 


3.33 


0.07 


2.13 


0.02 


-2.63 


0.06 


0.025 


0.10 


0.1210 




6668 


6:41:24.02 


9:26:53.0 


17.61 


0.01 


3.08 


0.03 


1.83 


0.01 


-2.76 


0.03 


0.008 


0.03 


0.0082 




6805 


6:41:25.93 


9:30:26.6 


17.89 


0.01 


3.45 


0.04 


2.11 


0.03 


-2.58 


0.05 


0.014 


0.05 


2.3883 


R- 


7037 


6:41:28.61 


9:53:58.5 


17.14 


0.01 


3.23 


0.02 


1.94 


0.02 


-2.64 


0.07 


0.013 


0.05 


12.7790 




7075 


6:41:29.16 


9:28:48.6 


19.39 


0.03 






2.37 


0.08 


-1.71 


0.10 


0.029 


0.11 


0.0000 




10865 


6:42:03.92 


9:49:02.2 


18.94 


0.02 


3.49 


0.13 


2.29 


0.03 


-2.64 


0.08 


0.037 


0.14 


45.2014 





